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This  document  describes  a  standard  procedure  for  empirically  de¬ 
termining  canopy  sun-glint  signatures,  a  standard  display  format 
for  the  sun-glint  signatures,  and  a  procedure  for  analytically 
deriving  canopy  sun-glint  signatures. 
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1.  Background 

Sun  glinting  from  the  canopy  is  a  major  cue  in  the  visual 
detection  of  aircraft.  An  aircraft  canopy  reflects  the  sun — 
produces  sun  glints — at  azimuth  and  elevation  angles  which  are 
a  function  of  the  geometry  of  the  canopy  and  the  orientation 
of  the  aircraft  with  respect  to  the  sun.  As  the  aircraft  ro¬ 
tates,  as  shown  in  Figure  1,  the  sun  glint  also  moves;  and,  as 
the  orientation  of  the  aircraft  changes  with  respect  to  the 
sun,  more  or  different  canopy  panels  can  become  the  reflecting 
surfaces.  At  any  one  instant,  canopy  sun  glints  may  be  visi¬ 
ble  at  more  than  one  set  of  azimuth  and  elevation  coordinates 
with  respect  to  the  aircraft. 


180° 


Figure  1.  Mechanism  of  Canopy  Sun  Glint 

To  reduce  the  chance  of  visual  detection,  aircraft  cano¬ 
pies  should  be  designed  such  that  no  sun  glints  fall  within  the 
band  of  coordinates  which  are  probable  positions  for  hostile 
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observers.  This  document  helps  to  meet  this  design  objective 
by  presenting: 

•  A  standard  experimental  procedure  for  determining 
canopy  sun-glint  signatures 

•  A  standard  format  for  documenting  canopy  sun-glint 
signatures 

•  An  analytical  procedure  for  determining  canopy  sun- 
glint  signatures 


2.  Experimental  Equipment 


Determination  of  canopy  sun-glint  signatures  is  performed 
using  a  scale  model  of  the  aircraft  being  studied,  a  simulated 
sun  source,  and  provisions  for  observing  and  recording  glint 
patterns.  A  typical  test  facility,  shown  in  Figure  2,  contains 
these  necessary  elements: 

•  Aircraft  model  •  Viewing-screen  support 

•  Model-mounting  fixture  •  Sun  simulator 

•  Viewing  screen  •  Sun-simulator  support 


Sun  Simulator 


Viewing  Screen 


Sun-Simulator 

Support 


Viewing-Screen  Model-Mounting  Aircraft  Model 

Support  Fixture 


Figure  2.  Test  Setup  for  Determining  Canopy  Sun-Glint  Signatures 
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AIRCRAFT  MODEL 


A  1/32-scale  model  is  convenient  for  these  tests.  If 
another  scale  is  selected,  the  diameter  of  the  viewing  screen 
must  be  changed  proportionally.  Surfaces  other  than  those 
being  studied  should  be  coated  with  flat  black  paint  to  elim¬ 
inate  extraneous  glints.  Since  this  procedure  for  determining 
sun-glint  signatures  does  not  measure  absolute  glint  intensi¬ 
ties,  only  position  and  coverage,  glint  intensity  may  be  en¬ 
hanced,  when  practical,  to  improve  visibility.  For  example, 
painting  the  interior  surfaces  of  transparent  canopies  with 
glossy  paint  increases  glint  intensity  without  changing  their 
size  or  shape. 


MODEL-MOUNTING  FIXTURE 


Any  convenient  mounting  fixture  with  a  head  which  rotates 
in  azimuth  can  be  used.  The  rotating  head  should  be  calibrated 
and  marked  in  degrees.  A  vertical-adjustment  capability  is 
required  in  either  the  model-mount  fixture  or  the  viewing- 
screen  support. 


VIEWING  SCREEN 


The  screen  must  be  translucent  to  permit  viewing  the 
projected  sun  glints  from  outside  the  cylinder.  Since  glint 
patterns  are  manually  marked  directly  on  the  exterior  surface 
of  the  screen,  the  screen  must  have  the  strength  and  the  abil¬ 
ity  to  accept  marking.  The  screen  must  be  replaceable  since  a 
new  screen  is  used  for  each  sun  elevation  angle  studied. 
Drawing  Mylar  has  been  used  successfully  as  a  screen  material. 
A  screen  diameter  of  48  inches  (1.22  meters)  is  used  for  1/32- 
scale  models. 


VIEWING-SCREEN  SUPPORT 


The  viewing-screen  support  must  be  sufficiently  rigid  to 
prevent  movement  with  respect  to  the  model  as  a  result  of  mark¬ 
ing  the  viewing  screen  while  outlining  reflected  sun  glints. 

The  screen  support  must  also  permit  installation  and  removal 
of  the  viewing  screens. 


SUN  SIMULATOR 


The  sun  simulator  is  an  incandescent  light  source,  aper- 
tured  by  a  long,  internally  absorbing  tube  to  provide  a  narrow 
beam  of  relatively  collimated  light.  The  beam  dispersion  must 
be  sufficient  to  cover  the  model. 
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SUN-SIMULATOR  SUPPORT 

The  sun-simulator  support  must  permit  movement  of  the  sun 
simulator  through  an  angle  in  the  vertical  plane  sufficient  to 
encompass  the  range  of  sun-elevation  angles  to  be  tested.  It 
must  be  possible  to  clamp  the  sun  simulator  at  desired  test  ele¬ 
vation  angles.  A  fixed  radius  from  the  model  is  not  required. 


3.  Experimental  Procedure 

The  model  is  mounted  on  its  mounting  fixture  with  the 
center  of  gravity  of  the  aircraft  on  the  nominal  rotation  axis 
of  the  mounting  fixture.  The  model  is  raised  or  lowered  until 
the  centerline  of  the  sun-simulator  axis  at  an  elevation  angle 
of  0  degrees  passes  through  the  nominal  eye  level  of  a  pilot 
seated  within  the  aircraft.  This  position  is  marked  as  the 
0-degree  line  on  the  viewing  screen.  The  sun  simulator  is  next 
positioned  at  the  desired  test  elevation  angle.  The  aircraft 
is  positioned  with  its  nose  pointed  at  the  sun  simulator  and 
its  waterline  parallel  to  the  viewing-screen  support  base.  The 
room  lights  are  extinguished  to  improve  visibility,  and  all 
glints  visible  on  the  viewing  screen  are  outlined  with  a  mark¬ 
ing  pen.  The  model  is  then  rotated  to  the  right  20  degrees  and 
the  glint-outlining  process  repeated.  The  sequence  of  rotating 
the  model  20  degrees  and  marking  sun  glints  is  repeated  until 
the  model  has  completed  360  degrees  of  rotation.  For  aircraft 
which  possess  left-right  symmetry,  0-  to  180-degree  rotation  of 
the  model  is  sufficient.  For  other  sun-elevation  angles,  the 
screen  is  replaced,  the  sun-elevation  angle  is  changed,  and  the 
procedure  is  repeated.  At  low  sun-elevation  angles,  a  gap  must 
be  left  in  the  viewing  screen  to  permit  illumination  of  the 
model.  It  may  be  necessary  to  extrapolate  for  glint  signatures 
reflected  at  the  0-degree  azimuth  in  such  cases. 


4.  Display  Format 


Figure  3  is  a  top  view  of  the  test  setup  shown  in  Figure 
2.  As  indicated,  the  sun  is  fixed  and  always  located  along  the 
0-degree  azimuth.  Any  observer  position  can  then  be  defined 
uniquely  in  terms  of  an  azimuth  and  elevation  measured  relative 
to  the  aircraft  rotational  axis  and  referenced  to  the  sun  posi¬ 
tion.  Aircraft  headings  are  also  referenced  to  the  sun  direc¬ 
tion  and  are  measured  positive  to  the  right.  If  canopy  sun 
glints  are  marked  on  the  screen,  then  the  outlines  of  the  sun 
glints  on  the  screen  define  the  azimuths  and  elevations  where 
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Figure  3.  Top  View  of  Experimental  Test  Setup 

sun  glints  can  be  seen  for  each  aircraft  heading.  If  the  view¬ 
ing  screen  is  opened  at  the  180-degree  azimuth  and  flattened 
with  the  outside  of  the  screen  face  up,  there  is  a  panoramic 
display  of  the  sun  glints  produced  by  the  aircraft  for  the  set 
of  heading  angles  selected.  This  display  defines  all  possible 
observer  locations  where  sun  glint  may  be  seen.  Fy  specifying 
the  markings  on  the  flattened  viewing  screen,  a  standard  dis¬ 
play  format  for  test  data  can  be  obtained.  Figure  4  is  a 
standard  display  developed  in  just  this  way. 


Figure  4  presents  the  canopy  sun-glint  signature  of  the 
AH-1G  in  standard  format.  This  display  contains  the  elements 
which  shall  be  included  in  the  documentation  of  canopy  sun- 
glint  signatures: 

•  Aircraft  and  canopy  identification 

•  Sun-elevation  angle  and  aircraft  pitch  and  roll 
attitude 

•  Scale  ranging  from  0  to  -180  and  0  to  180  degrees  of 
azimuth  locations  from  the  aircraft,  referenced  to 
the  direction  of  the  sun 
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•  Analyst-specified  band  of  potential  observer  locations 
with  respect  to  the  aircraft 

•  Canopy  sun-glint  signatures,  together  with  the  air¬ 
craft  heading  with  respect  to  the  sun  which  produced 
them 

•  Probability  of  sun  glint  occurring  in  the  display  area 
bounded  by  the  observer  locations  and  the  0  to  -180 
and  0  to  180-degree  azimuth  boundaries 


Figure  4.  Standard  Format  for  Presenting  Canopy  Sun-Glint  Signature 
SUN-ELEVATION  ANGLE  AND  AIRCRAFT  PITCH  ATTITUDE 


The  sun  elevation  angle  shall  be  presented  above  and  on 
the  left  side  of  the  display  as  shown  in  Figure  4  (A) .  Air¬ 
craft  pitch  attitude  and  roll  attitude  shall  be  presented 
directly  below  the  sun  elevation  angle. 


AZIMUTH 

The  horizontal  space  devoted  to  the  display  shall  be 
divided  into  twelve  equal  segments  as  shown  in  Figure  4  (B) 

and  the  azimuth  position  of  the  sun  glints  with  respect  to  the 
sun  marked  in  30-degree  increments  at  the  bottom  of  the  dis¬ 
play.  The  display  is  viewed  as  an  opened  right  circular 
cylinder  with  0  degrees  being  in  the  direction  of  the  sun  and 
the  angles  increasing  positively  to  the  right  and  negatively  to 
the  left  as  seen  from  the  inside  of  the  cylinder.  The  cylinder 
axis  is  vertical 

POTENTIAL  LOCATIONS  OF  OBSERVERS 

The  observer-elevation  angles  shall  be  shown  as  in  Figure 
4  (C )  .  The  minimum  observation  band  shall  be  ±10  degrees  from 
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the  aircraft  horizon.  The  eye  level  of  the  aircraft  pilot  is 
nominally  taken  as  0  degrees. 


AIRCRAFT  HEADING 


With  two  exceptions,  each  canopy  sun  glint  which  falls 
within  the  display  angular  boundaries  as  the  aircraft  is  ro¬ 
tated  on  its  yaw  axis  in  20-degree  increments  shall  be  shown 
as  in  Figure  4  (D) .  If  the  canopy  sun-glint  signature  is  so 
extensive  that  rotating  the  aircraft  causes  considerable  over¬ 
lapping,  a  single  heading  may  be  selected  for  display.  (The 
OH-6  data  of  this  report  is  an  example  of  this  exception.) 


If  the  canopy  signature  changes  little  with  aircraft 
rotation,  other  than  a  translation,  data  may  be  shown  in  incre¬ 
ments  of  40  degrees.  (The  OH-58  data  of  this  report  is  an  ex¬ 
ample  of  this  exception.) 

Aircraft  heading  shall  be  given  as  0  degrees  with  the 
nose  pointed  at  the  sun.  Rotation  shall  be  positive  to  the 
right.  For  right-left  symmetric  canopies,  only  the  sun  glints 
generated  for  0  to  180  degrees  of  aircraft  heading  rotation 
need  be  displayed.  Further  rotation  would  only  yield  the 
mirror  image  of  the  glint  signatures  generated  in  the  first 
0  to  180  degrees  of  rotation. 


PROBABILITY  OF  SUN  GLINT 


For  each  display,  the  probability  of  a  glint,  Figure  4 
(E) ,  existing  within  the  observer  elevation  band  shall  be  com¬ 
puted.  Probability  of  glint  is  defined  as  the  fraction  of  the 
panorama  area  between  0  and  -180  degrees  and  0  and  180  degrees 
with  respect  to  the  sun  and  between  the  upper  and  lower  observ¬ 
er  elevation  angle  limits  which  is  swept  by  the  reflected 
canopy  sun  glint  as  the  aircraft  is  rotated  360  degrees  in 
heading.  A  swept  area  is  counted  only  once.  For  example,  if 
the  area  between  60  degrees  and  90  degrees  with  respect  to  the 
sun  and  0  degrees  and  10  degrees  of  observer  elevation  angle 
is  the  position  of  a  glint  with  aircraft  heading  at  both  20 
degrees  and  340  degrees,  the  fact  that  there  is  an  overlap  is 
ignored  in  the  computation.  For  the  limits  shown  on  Figure  1, 
there  is  360  degrees  x  20  degrees  =  7,200  degrees2  total  area. 
The  sun  glint  sweeps  over  305  degrees  x  15  degrees  or  4,575 
degrees2.  The  probability  of  glint  is  therefore  4,575/7,200  = 
0.6354.  The  glint  probability  shall  be  displayed  below  and  to 
the  right  of  the  glint  panorama  as  shown  in  Figure  4  (E) . 
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5.  Experimentally  Derived  Canopy  Sun-Glint  Signatures 


The  following  four  figures  display  canopy  sun-glint  sig¬ 
natures  of  four  different  canopy  configurations  in  standard 
format.  Sun-elevation  angles  were  chosen  to  bracket  the  sun- 
elevation  angles  occurring  in  the  Fulda  Gap  area  of  the  Federal 
Republic  of  Germany.  Figure  5  is  the  canopy  sun-glint  signa¬ 
ture  of  the  AH-1G  with  the  standard  canopy.  Figure  6  is  the 
canopy  sun-glint  signature  of  the  AH-1G  equipped  with  a  canopy 
having  vertical  sides,  flat  top  panel  and  front  windshield,  and 
a  canopy  centerline  fence.  Figures  7  and  8  display  the  canopy 
sun-glint  signatures  of  the  OH-6  and  OH-58  respectively. 


6.  Analytical  Procedure 


BASIC  CALCULATIONS 

The  sun  glint  from  a  reflecting  surface  depends  upon  the 
reflection  angles  of  the  sun's  rays  from  each  boundary  point 
of  the  surface.  Determination  of  these  angles  depends  on  the 
angle  of  incidence  of  the  sun's  rays  and  the  normal  to  the 
reflecting  surface.  The  angle  of  incidence  is  given  by  the 
specified  sun  elevation  and  sun  azimuth,  while  the  normal  to 
the  surface  depends  on  the  type  of  surface  and,  in  general, 
the  coordinates  of  any  point.  For  example,  in  the  case  of  a 
geometric  plane,  or  flat  surface,  the  normal  to  the  surface 
has  the  same  direction  at  every  point.  However,  the  quadratic 
or  second-order  surface  has  a  normal  which  varies  in  direction 
from  point  to  point.  These  two  types  of  surfaces  were  consid¬ 
ered  to  be  representative  of  most  canopy  shapes.  Combinations 
of  these  two  types  would  generate  more  complex  shapes. 


Surface  normals  and  the  angles  of  incidence  must  be  re¬ 
ferred  to  a  common  axis  system  in  order  to  permit  the  calcula¬ 
tion  of  reflection  angles.  A  convenient  axis  system  is  one 
which  is  fixed  in  space  (see  Figure  9) .  The  positive  Z-axis 
points  toward  the  center  of  the  earth.  The  X-axis  is  chosen 
perpendicular  to  the  Z-axis  and  in  such  a  direction  as  to  cause 
the  sun  to  lie  in  the  X-Z  plane.  In  other  words,  a  sun  eleva¬ 
tion  angle  of  less  than  90  degrees  means  that  the  sun  has  a 
positive  X  and  negative  Z  direction.  It  will  also  have  zero 
azimuth  referred  to  this  axis  system.  Finally,  the  Y-axis  is 
chosen  so  as  to  complete  a  right-handed  coordinate  system. 
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Figure  5.  Canopy  Sun-Glint  Signature  of  AH-1G  With  Standard  Canopy 
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Figure  7.  Canopy  Sun-Glint  Signature  of  OH-6  With  Standard  Canopy 
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Figure  9.  Fixed-Coordinate  System  for  Angular  Measurements 


Since  the  reflecting  surface  is  normally  described  in 
terms  of  an  axis  system  fixed  to  the  aircraft,  a  transforma¬ 
tion  of  points  and  curves  from  one  axis  system  to  the  other  is 
needed.  This  transformation  of  body-fixed  axes  to  earth-fixed 
axes  consists  of  three  aircraft  rotations  (roll,  pitch,  and 
yaw) .  It  is  assumed  that  the  origins  of  both  axis  systems 
coincide. 

Next,  if  there  exist  one  or  more  fences  on  the  aircraft, 
the  effect  of  shadows  cast  by  these  fences  must  be  analyzed. 
This  is  done  by  projecting  the  boundary  points  from  both  the 
fences  and  reflecting  surface  along  the  sun's  rays  onto  a 
plane  perpendicular  to  these  rays.  The  intersection  of  the 
shadow  and  the  reflecting  surface  on  this  plane  will  determine 
a  new  set  of  boundary  points  for  the  reflecting  surface.  This 
new  set  of  points  will  represent  the  portion  of  the  surface 
remaining  outside  of  the  shadow. 

The  normal  to  the  surface  at  each  remaining  or  new  bound¬ 
ary  point  is  then  calculated  based  on  the  geometric  shape  of 
the  surface.  Generally,  it  is  assummed  that  a  great  deal  is 
known  about  the  reflecting  surface.  Either  an  exact  mathe¬ 
matical  representation  or  general  shape  should  be  known.  How¬ 
ever,  if  this  is  not  so,  the  boundary  points  may  be  curve- 
fitted  to  a  specified  type  of  surface  using  a  modified  least - 
squares  error  procedure.  This  surface  equation  then  enables 
the  rapid  determination  of  surface  normals. 

As  was  mentioned,  the  reflection  vectors  from  each  bound¬ 
ary  point  are  determined  using  the  surface  normal  and  angle  of 
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incidence  of  the  sun’s  rays.  Assuming  a  smooth  reflecting  sur¬ 
face,  the  reflection  vector  is  equal  to  the  sun  vector  with 
its  normal  component  to  the  surface  negated.  It  must  now  be 
determined  if  the  reflection  vector  intersects  any  fences  on 
its  way  to  the  observer.  This  is  done  by  extending  the  vector 
until  it  intersects  the  plane  of  the  fence.  If  this  inter¬ 
section  lies  within  the  boundary  points  of  the  fence,  then  the 
observer  will  not  view  the  reflection. 

Using  all  the  remaining  reflection  vectors,  a  sun  glint 
at  an  assumed  observer  distance  can  be  calculated.  This 
glint  is  referred  to  a  specific  reference  point.  Determina¬ 
tion  of  this  sun  glint  is  obtained  by  extending  each  reflec¬ 
tion  vector  until  it  intersects  a  cylinder  whose  radius  equals 
the  distance  to  the  observer.  By  measuring  the  elevation  and 
azimuth  angles  of  the  line-of-sight  vector  from  the  reference 
point  to  each  intersection,  a  sun  glint  is  formed. 


DIGITAL  COMPUTER  PROGRAM 


A  digital  computer  program  was  written  to  solve  in  detail 
the  above-described  glint  problem.  It  was  coded  in  FORTRAN  IV 
and  has  been  made  operational  on  the  IBM  360  system.  The  main 
(executive)  routine  controls  the  general  flow  of  the  program. 
It  reads  in  the  desired  helicopter  angles,  plot  information, 
observer  distance,  and  boundary  points  in  body  axes  of  both 
fences  and  reflecting  surfaces.  It  sets  up  necessary  DO  loops 
and  calculates  surface  normals.  Figure  10  is  a  flow  chart  for 
the  main  program. 

A  general  description  of  all  the  subroutines  used  in  the 
entire  program  follows: 

•  CFIT — Sets  up  boundary  points  in  desired  form  to  be 
curve  fitted. 

•  CURFIT — Sets  up  necessary  matrices  for  least-square 
curve  fitting  of  given  points.  Calculates  standard- 
deviation  error  for  calculated  fit. 

•  EQNSOL — Solves  a  system  of  linear  equations,  either 
homogeneous  or  inhomogeneous.  This  subroutine,  along 
with  CFIT  and  CURFIT,  curve  fits  boundary  points  to 
either  a  plane  or  a  quadratic  surface,  depending  on 
type  specified  by  input. 

•  CVBI — Uses  the  body-to-inertial-axes  transformation  to 
convert  coefficients  of  our  surfaces  in  body  axes  to 
coefficients  in  inertial  axes. 
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•  SHADOW — Determines  the  intersection  of  a  shadow  from  a 
fence  with  reflection  surface.  Defines  new  reflection 
surface  boundary  points  based  on  portion  of  surface 
outside  of  shadow. 

•  PLANPT — Calculates  the  intersection  of  a  line  with  a 
plane.  The  inputs  are  a  point  in  space  on  the  line 
and  direction  numbers  for  the  line  along  the  coeffi¬ 
cients  of  the  plane. 

•  INTERC — Finds  intersections,  if  any,  of  a  line  given 
by  two  points  in  a  plane  with  set  of  lines  given  by  a 
set  of  connected  points  in  same  plane. 

•  CHECK 3 — Checks  to  see  whether  a  particular  set  of 

X  and  Y  values  of  a  point  lie  outside  of  a  given  range. 

•  SORNOT — Determines  whether  a  particular  boundary  point 
of  a  fence  is  between  sun  and  reflection  surface  or 
behind  it.  It  is  also  used  in  regard  to  finding  inter¬ 
section  of  reflection  vector  with  cylinder  of  radius 
equal  to  observer  distance.  In  this  latter  case  it  is 
used  to  choose  between'  two  possible  solutions  by  deter¬ 
mining  which  solution  lies  in  positive  direction  of 
reflection  vector. 

•  TRANSP — Contains  transformation  that  takes  inertial 
axes  into  set  of  axes  where  the  Z-axis  lies  in  the 
same  direction  as  the  sun  vector.  If  a  set  of  points 
lies  in  a  plane  perpendicular  to  the  sun,  this  sub¬ 
routine  converts  this  set  of  three-dimensional  inertial 
coordinate  of  points  to  a  set  of  two-dimensional  coordi¬ 
nates  of  points  in  this  new  axis  system.  It  also  pro¬ 
vides  inverse  for  transforming  points  back  to  inertial 
axes. 

•  CLOSE — Determines  which  point  in  set  of  points  is 
closest  to  a  given  point. 

•  CHECKl — Determines  whether  a  given  point  is  already 
contained  in  a  set  of  points. 

•  CHECK2 — Similar  in  function  to  CHECKl  except  that  the 
given  point  must  first  be  transposed  to  same  coordi¬ 
nate  system  that  set  of  points  is  in. 


•  CIRCLE — Determines  whether  a  given  point  is  enclosed 
by  a  given  set  of  points  if  each  point  in  the  set  is 
connected  by  a  straight  line. 
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•  INTERF — Calculates  whether  reflection  from  any  par¬ 
ticular  boundary  point  intersects  any  fence. 

•  TRANQ — Projects  points  on  quadratic  surface  to  points 
in  plane,  where  the  plane  is  a  best-fit  plane  through 

a  set  of  quadratic  surface  boundary  points.  Projection 
is  along  sun  rays.  This  subroutine  will  also  perform 
reverse  projection  assuming  that  nearest  intersection 
of  quadratic  surface  preserves  one-to-one  mapping. 

•  QUADPT — Calculates  intersection  or  intersections  of  a 
line  with  a  quadratic  surface.  Inputs  to  the  sub¬ 
routine  are  a  point  on  a  line,  direction  numbers  of 
the  line,  and  coefficients  of  the  quadratic  surface. 

•  PLOT — An  on-line  digital  plot  of  glints  from  a  re¬ 
flecting  surface  for  a  range  of  helicopter  headings 
for  a  given  sun  elevation  and  helicopter  pitch  and 
roll  attitudes. 

•  PROBL — Calculates  swept-out  area  between  glints  on  an 
elevation-versus-azimuth  basis.  Minimum  and  maximum 
elevation  angles  are  used  with  straight  lines  drawn 
between  corresponding  values  for  different  glints. 

Area  is  calculated  using  triangular  area  law  of 
points. 

•  ADDPTS — Will,  if  desired,  supply  additional  boundary 
points  along  curved  portions  of  quadratic  surfaces. 

A  detailed  description  of  all  input  variables  and  their 
format  is  contained  in  Appendix  A.  A  complete  FORTRAN  IV 
listing  of  the  program  appears  in  Appendix  B.  Basically  the 
inputs  are  desired  aircraft  attitudes,  sun  elevations,  ob¬ 
server  elevations,  and  boundary  points.  An  assumption  was 
made  that  all  fences  could  be  handled  as  combinations  of  flat 
surfaces  only.  For  the  purpose  of  calculating  shadows,  it  was 
felt  that  this  was  sufficient.  As  mentioned  before,  reflect¬ 
ing  surfaces  were  assummed  to  be  combinations  of  flat  and 
quadratic  shapes.  It  was  also  necessary,  in  order  to  save 
some  sorting  and  execution  time,  to  require  that  the  boundary 
points  be  simply  connected.  That  is,  point  1  is  connected  to 
point  2  and  point  2  is  connected  to  point  3,  and  so  on,  with 
the  last  point  connected  to  point  1.  No  other  combination  is 
allowed. 

In  calculating  the  intersection  of  a  shadow  with  a  re¬ 
flecting  surface,  the  assumption  was  made  that  the  points 
could  be  connected  by  straight  lines.  An  input  to  the  program 
will  cause  the  generation  of  additional  boundary  points  for 
curved  portions,  if  necessary,  thus  minimizing  the  error  re¬ 
sulting  from  this  straight-line  assumption. 
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An  example  of  a  complete  set  of  inputs  is  tabulated  in 
Figure  11.  These  inputs  are  for  the  AH-1G  with  straight  sides. 
The  configuration  is  symmetric  with  respect  to  the  X-Z  plane, 
so  only  the  top,  front,  and  left  side  are  needed  to  generate 
the  glint  signature.  Also,  only  180  degrees  of  aircraft  rota¬ 
tion  is  needed. 

Output  to  the  program  is  in  the  form  of  tables  and  on-line 
plots,  as  shown  in  Figures  12  through  16.  The  tables  consist 
of  all  necessary  input  information  along  with  reflection  vector 
and  glint  angles.  Figure  12  is  a  tabulation  of  boundary  points 
for  the  center  fence  and  the  resulting  curve-fit.  Figure  13  is 
a  tabulation  of  boundary  points  for  the  first  reflection  surface 
and  the  resulting  curve-fit.  Figure  14  is  a  tabulation  of  the 
output  reflection  vector  and  glint  angles  as  a  function  of  air¬ 
craft  azimuth  for  constant  sun  elevation  and  aircraft  pitch  and 
roll.  This  table  is  for  reflection  surface  number  one.  The  in¬ 
formation  contained  in  Figures  13  and  14  is  repeated  for  each 
reflection  surface.  Also,  if  on-line  plots  are  requested.  Fig¬ 
ure  15  shows  the  plot  for  the  first  reflecting  surface  corres¬ 
ponding  to  the  data  of  Figure  14.  There  will  be  an  on-line 
plot  for  each  surface  at  each  sun  elevation  with  the  numbered 
points  corresponding  to  the  different  input  aircraft  yaw  azi¬ 
muths,  as  shown  in  Figure  15.  In  addition  to  a  plot  for  each 
separate  surface,  there  will  be  a  composite  plot  of  all  surfaces 
for  each  desired  sun  elevation,  as  is  shown  in  Figure  16. 


VERIFICATION  OF  ANALYTICAL  PROCEDURE 


The  first  step  in  verifying  the  above-described  computer 
program  was  to  input  simple  shapes  with  known  solutions.  Once 
this  was  done  successfully,  the  next  step  was  to  input  canopy 
configurations  for  which  experimental  data  was  available. 
Results  for  three  configurations  were  compared.  The  three 
configurations  were  the  AH-1G  with  vertical  sides,  the  OH-6, 
and  the  OH-58. 

Figure  17  is  the  glint  signature  for  the  AH-1G  with  verti¬ 
cal  sides  as  shown  in  Figure  6  with  the  computer  results  spotted 
on.  As  can  be  seen,  the  comparison  is  good.  This  is  to  be 
expected  since  all  surfaces  are  planes  and  thus  have  constant 
normal  directions.  Some  discrepancies  appear  in  the  range  of 
azimuths  corresponding  to  glints  resulting  from  side  window 
reflections  at  zero  sun  elevation.  There  are  at  least  three 
possible  causes.  The  plastic  canopy  used  in  generating  the 
experimental  data  was  not  smooth,  causing  some  dispersion.  The 
light  source  used  in  the  experiment  had  some  associated  disper¬ 
sion.  The  third  error  could  be  in  the  computer  input  data,  in 
that  the  dimensions  were  not  generated  from  a  detailed  drawing 
but  calculated  from  a  partially  scaled  drawing. 


20 


i 


Input 

Card 

Type 


wbwjwk.^'w -"vmx?***' 


Input  Data 


& 

i 


1 

•180 , 

180. 

-100. 

100. 

2 

0.00 

0.0 

0.0 

60. 

•10 

10. 

3 

•  0.0 

0.0 

4 

6.0 

5. 

2. 

6 

10 

0,00 

10.0 

20,0 

30.0 

00,0 

7 

50. 

60, 

70. 

80, 

90. 

7 

100. 

120, 

1O0. 

160. 

8 

1 

9 

.0 

10 

13.170 

0. 

1. 

11. 

0. 

-1.88 

10 

o.o 

0. 

•1.88 

0. 

o. 

1. 

11 

0 

0 

1 

1  0 

12 

1 

0 

0 

0 

13 

12.79 

•1.5 

1. 

12.79 

1.5 

1. 

13 

11. 

1.5 

•1,38 

11. 

-1.5 

•1.38 

12 

1 

0 

0 

0 

13 

11. 

•1.5 

•1.38 

11. 

1.5 

•1,38 

13 

«. 

1.5 

•1,38 

0. 

•1.5 

•1,38 

12 

1 

3 

0 

0 

13 

12.79 

•1.5 

1. 

11. 

•1.5 

-1.38 

13 

11. 

•1.5 

1.0 

12 

1 

0 

0 

0 

13 

11. 

•1.5 

•1.38 

9. 

•1 .5 

-1,38 

13 

«. 

•1.5 

1. 

11. 

•1.5 

1. 

Figure  11.  Computer  Printout 

of  Inputs  for 

the  AH-1G 

With  Flat- 

Side  Canopy 

and  Centerline  Fence 


21 


Figure  14.  Computer  Printout  of  Output  Reflection  and  Glint  Angles 
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Figure  17.  Computer  Results  Superimposed  on  Sun-Glint  Signature 
of  AH-1G  With  Flat-Side  Canopy  and  Centerline  Fence 


26 


Figure  18  is  the  superimposed  computer  results  for  the 
OH-6  at  specified  aircraft  azimuths.  Here,  as  can  be  seen, 
the  digital  results  are  not  nearly  as  good  as  the  AH-1G  with 
vertical  sides.  Possible  causes  are: 

•  The  plastic  model  may  not  have  been  completely  repre¬ 
sentative  of  the  OH-6.  Slight  changes  in  certain  di¬ 
mensions  can  cause  apparently  large  changes  in  results. 
For  instance,  a  change  of  0.05  inch  in  the  model  cut¬ 
out  of  the  front  windshield  could  cause  a  20-degree 
shift  in  azimuth  of  one  of  the  glint  boundaries. 

•  The  exact  location  of  the  reference  point  used  in  the 
generation  of  the  experimental  data  was  not  known. 

•  The  mathematical  representation  for  the  reflecting 
surfaces  was  obtained  from  three-view  drawings  with 
some  necessary  simplification. 

However,  the  general  location  and  shape  of  the  glints  were  in 
agreement.  It  must  also  be  remembered  that  a  1-degree  error 
in  calculation  of  surface  normal  will  result  in  a  2-degree 
error  in  the  position  of  the  reflection.  For  rapidly  changing 
normals,  this  multiplying  factor  becomes  significant. 

Figure  19  is  the  glint  signature  for  the  OH-58  with  com¬ 
puter  results  plotted  on  top.  Of  the  three  configurations, 
the  mathematical  expression  of  the  windshield  of  the  OH-58  was 
the  most  difficult  to  determine  and  thus  the  most  questionable. 
Also, the  question  of  model  exactness  must  be  raised  again. 

The  windshield  is  not  a  simple  shape,  and  without  detailed 
cross-sectional  drawings,  the  representation  of  the  windshield 
as  a  combination  of  simple  shapes  was  difficult.  The  refer¬ 
ence  point  again  was  not  known.  However,  as  with  the  OH-6, 
the  results  show  general  agreement. 

The  conclusion  from  these  comparisons  is  that  if  one  is 
looking  for  relative  merits  or  problems  of  a  particular  canopy 
design,  then  curve-fitting  of  relatively  crude  inputs  will 
yield  adequate  results.  However,  if  absolute  answers  in  terms 
of  precise  location  of  each  glint  corresponding  to  each  air¬ 
craft  and  sun  location  are  desired,  then  a  very  good  grasp  is 
needed  on  the  mathematical  expression  of  the  canopy  design. 

This  last  statement  is  not  as  confining  as  it  may  seem,  for  if 
a  person  is  designing  a  canopy,  he  will  normally  know  the 
precise  shape  of  his  canopy. 
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Figure  19.  Computer  Results  Superimposed  on  Sun-Glint  Signature 
of  OH-58  With  Standard  Canopy 
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Appendix  A 

Input  and  Output  of  Digital-Computer  Program 


INPUT 


All  card  types  with  (*)  are  optional,  depending  on  the  type 
and  number  of  inputs. 


Card  Type  1 — Format  4E10.0 


XSCI1IN 

Lower  limit  for  glint  azimuths  for  on-line 
plots 

XSCMAX 

Upper  limit  for  glint  azimuths  for  on-line 
plots 

YSCMIN 

Lower  limit  for  observer  angles  for  on-line 
plots 

YSCMAX 

Upper  limit  for  observer  angles  for  on-line 
plots 

Card 

Type  2- 

-Format  6E10.0 

XREF 

X  distance  along  X-body  axis  to  reference 
point,  reference  point  being  point  from  which 
glint  angles  are  measured 

YREF 

Y  distance  to  reference  point 

ZREF 

Z  distance  to  reference  point 

DISTG 

Distance  in  X-Y  inertia  plane  to  assumed 
observer 

OBSMIN 

Lower  limit  on  observer  elevation  angle  used 
in  calculating  probability  function 

OBSMAX 

Upper  limit  on  observer  elevation  for  prob¬ 
ability  calculation 

Card 

Type  3- 

-Format  2E10.0 

THETA 

Pitch  attitude  of  aircraft  (positive  up) 

PHI 

Roll  attitude  of  aircraft  (positive  right 

wing  down) 
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Card  Type  4 — Format  15,5X,5E10.0 

NGAM  Number  of  sun-elevation  angles  for  which  data 

is  desired 

GAM  Array  containing  desired  sun  elevations 

Card  Type  5 — Format  5E10.0 

GAM  If  NGAM  >  5,  additional  values  of  sun 

elevations 

Card  Type  6 — Format  15 f  5<,  5E10.0 

NPSI  Number  of  aircraft  azimuths  for  which  data 

is  derired 

PS  Array  containing  desired  aircraft  azimuths 


; 


i 

i 

* 


i 

i 

i 


>■ 


Card  Type  7 — Format  5E10.0 

PS  If  NPSI  >  5,  additional  values  of  aircraft 

azimuths 

Card  Type  8 — Format  15 

I 

NFENCE  Number  of  fences  for  which  boundary  point 
will  be  read  in 

Card  Type  9* — Format  15 

NFPTS  Array  containing  number  of  boundary  points 

for  each  fence 

Card  Type  10* — Format  6E10.0 

I 

FPTS  Array  containing  boundary  points  for  each 

fence.  X,Y.Z  locations  in  body  axis  are  read 
in  with  two  points  on  a  card.  This  card  type 
is  continued  to  be  read  until  all  points  for  a 
particular  fence  are  read  in.  Then  another 
type-9  card  is  read,  and  so  on,  until  all 
fences  are  read  in.  Card  types  9  and  10  are 
used  only  if  NFENCE  >  0. 

Card  Type  11 — Format  515 

NPANEL  Number  of  reflecting  surfaces  to  be  read  in 

IANGL  Determines  form  of  glint  azimuths  in  print 

out: 


IANGL  =0  ±180  degrees 

IANGL  =1  0  degrees  to  360  degrees 
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IPLOT  Are  on-line  plots  desired? 

I PLOT  =0  No 

IPLOT  =  1  Yes 

IPROB  Is  probability  function  to  be  calculated? 

IPROB  =0  No 

IPROB  =  1  Yes 

NDISP  Is  sun  dispersion  of  0.5  degrees  to  be  used? 

NDISP  =  0  No,  parallel  rays 

NDISP  =  1  Yes 


Card  Type  12 — Format  515 


ITYPE 


NPTS 

ICOEFS 


NBETW 


Type  of  reflecting  surface: 


ITYPE  =  1  Flat  surface 

ITYPE  =  2  Quadratic  surface 

Number  of  boundary  points  for  reflecting 
surface 


Determines  whether  coefficients  for  reflecting 
surface  will  be  read  in  or  curve-fitted  from 
boundary  points.  Also  can  be  used  to  specify 
shape  of  fit: 


ICOEFS  =  0  Curve  fit  data  using 
built-in  method  of 
selecting  terms 

ICOEFS  =  1  Read  in  coefficients 
ICOEFS  =  2  Particular  terms  to  be 

fitted  will  be  specified 


Are  additional  boundary  points  along  curved 
surfaces  desired? 


NBETW  =0  No 

NBETW  =  N  Yes .  N  =  Number  of  addi¬ 

tional  points  between  each 
pair  of  points  read  in 

INTERN  Are  internal  cockpit  reflections  permitted? 


INTERN  =  0  Yes 

INTERN  =1  NO 
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Card  Type  13 — Format  6E10.0 

XPTSB  Array  containing  boundary  points  of  reflecting 

surface  in  body  axis.  X,Y,Z  components  of 
each  point  are  read  in  with  two  points  to  an 
input  card.  This  card  type  is  continued  to  be 
read  in  until  all  points  are  used. 


I 


l. 

ft 


Card  Type  14* — Format  1015 

NCOEFS  Number  of  coefficients  in  desired  curve  fit 

IFIT  Array  containing  location  of  each  desired  term 

in  basic  equation.  This  card  type  used  only 
if  ICOEFS  =  2  and  ITYPE  =  2.  The  following  is 

an  example  of  how  this  is  inputed.  The  basic 

quadratic  surface  equation  as  assumed  by  the 
program  is 

AX  Z+BY  2+CZ  2+DXY+EXZ+FYZ+GX+HY+KZ+L  =  0 

If  it  is  desired  to  curve  fit  boundary  points 
to  the  equation 

AX2+CZ2+HY+KZ+L  =  O 

Then  the  inputs  would  be  NCOEFS  =  4, 

IFIT  (1)  =  1,  IFIT  ( 2 )  =  3,  IFIT  ( 3)  =  8, 

IFIT  ( 4)  =  9. 

Constant  term  is  not  counted.  Also, 

IFIT  must  contain  integers  in  ascending 
magnitudes  only. 

Card  Type  15* — Format  6D10.0 

COEFB  Array  containing  coefficient  of  reflecting 

surface  equation 

For  ITYPE  =  1,  COEFB ( 1 ) X  +  COEFB(2)Y 

+  COEFB (3) Z  +  COEFB (4)  =  0 

For  ITYPE  =  2,  COEFB ( 1 ) X 2  +  COEFB ( 2) Y 2 
+  COEFB ( 3) Z  2  +  COEFB (4)XY 
+  COEFB ( 5 ) XZ  +  COEFB ( 6) YZ 
+  COEFB ( 7) X  +  COEFB ( 8) Y 
+  COEFB ( 9) Z  +  COEFB (10)  =  0 

This  card  type  is  used  only  if  ICOEFS  =  1. 

Two  cards  of  this  type  are  required  to  input 
the  10  coefficients.  Zeros  must  be  included 
for  undesired  coefficients. 
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OUTPUT 


The  output  to  the  program  is  simply  a  tabulation  and,  if 
desired,  an  on-line  plot  of  the  elevation  and  azimuth  angles  of 
each  glint  boundary  point  associated  with  each  reflection  sur¬ 
face  boundary  point.  Reflection  boundary  points  which  are 
read  in  but  lie  in  a  shadow  of  a  fence  will  not  appear  in  out¬ 
put,  and  points  whose  reflection  vector  intersects  a  fence  will 
have  'reflection  vector  angles  and  glint  angles  printed  as  a 
row  of  asterisks. 

The  elevation  and  azimuth  angles  are  calculated  for  each 
desired  sun  elevation  and  aircraft  yaw  angle,  and  are  referred 
to  a  particular  reference  point.  An  additional  output,  if  re¬ 
quested,  is  a  simplified  estimate  of  the  area  covered  by  the 
sun  glints  as  the  aircraft  rotated  through  a  range  of  azimuths. 
This  is  put  in  terms  of  a  ratio  to  total  area  and  is  termed 
the  "probability  function" . 

\ 
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Appendix  B 
Program  Listing 

SUGGESTIONS  FOR  PROGRAI1  USERS 


The  analytical  procedure  given  in  this  report  is  sensitive 
to  errors  in  canopy  geometry.  It  is  suggested  that  the  user 
take  advantage  of  defined  fuselage  lines  from  engineering  de¬ 
sign  data  when  available  for  the  aircraft  canopy  under  consid¬ 
eration.  For  flat-panel  canopies,  this  degree  of  precision  is 
not  required. 

The  analytical  procedure  is  designed  to  operate  using 
simply  connected  boundary  points  of  the  individual  canopy 
panels.  For  curved  canopies,  it  is  suggested  that  an  initial 
run  be  made  with  boundary  points  and  nonboundary  points  in¬ 
cluded  in  order  to  obtain  a  good  equation  for  the  surface  of 
the  canopy.  For  this  purpose  request  a  0-degree  sun  elevation 
and  a  0-degree  aircraft  heading  and  check  the  output  to  see  if 
the  reflection  vectors  are  in  the  expected  direction.  For 
example,  if  the  canopy  is  parabolic  so  that  a  normal  to  the 
highest  point  on  the  canopy  is  perpendicular  to  the  aircraft 
waterline,  then  the  reflection  vector  should  be  at  nominally 
180  degrees  azimuth  and  a  low  elevation — 0  to  10  degrees. 
Similarly,  if  a  normal  to  the  outermost  point  on  the  canopy  in 
the  Y  direction  is  perpendicular  to  the  X-Z  plane,  the  reflec¬ 
tion  vector  should  be  at  nominally  180  degrees  azimuth  and  0 
degrees  elevation.  By  continuing  this  procedure  for  other 
points  on  the  canopy  where  the  normal  to  the  canopy  is  known, 
the  user  can  assure  himself  that  the  equation  for  the  canopy 
is  reasonable.  As  a  starting  point,  specifying  only  equation 
terms  2,  3,  7,  and  9  on  card  type  14  has  been  found  useful  for 
parabolic  canopies.  Once  the  canopy  equation  has  been  defined, 
the  user  can  remove  the  nonboundary  points  from  the  data  set 
and  specify  the  canopy  equation  terms  for  further  runs.  The 
added  runs  would  encompass  all  sun  positions  and  aircraft 
positions  desired. 

There  is  a  note  of  warning  when  running  with  fences.  Bound¬ 
ary  points  which  were  inputted  may  not  appear  in  the  output; 
boundary  points  which  were  not  inputted  may  appear  in  the  out¬ 
put.  This  is  due  to  the  fact  that  a  new  set  of  boundary  points 
is  generated  whenever  the  shadow  of  the  fence  intersects  the 
reflective  surface.  Additional  boundary  points  may  also  appear 
in  the  output  when  the  shadow  does  not  intersect  the  reflective 
surface  if  the  reflective  surface  is  curved  (I TYPE  =  2) .  This 
occurs  because  the  intersection  of  the  sun's  vector  with  the 
curved  surface,  under  a  given  condition,  may  have  two  solutions. 
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(DIMENSION  FPT3(5,30,3),  NRPTS(5)  ,RPSI  (5,30,3) , 

1  PPSI(30#3)#NFPTS(  5 ) » GAM (4 )  » PS ( 16) 

(DIMENSION  COEFI(10)#OIRC08(3,3)#VNORM(3),SUN(3),REFLTN(3) 
DIMENSION  XPT8R ( 30  #  3 ) 

data  tol/i,oe«o«/ 

REAL*8  SIG,C0EP'B(10)  ,COEFD(10) 

REAl*8  CoEFHQ(IO) #C0EFBF(5,4) 

DIMENSION  COEFSUO)  ,COtFIPCtO) 
COMMON/CTR02/OIRCOS#ypT8R,XMIN(3)*XMAX(3)#NPTS 
REAL'S  OlRCOS,COEFS,CaEFI(COEFIF,8tlN,REFLTN,DI3C,X,  Y,Z»X1 , Y1 ,Z1 , 
1  Xt2, VI?,ZI2,0Q2,DATAn2,OAB8,D8QRT 
COMMON/E ENCE/C0EF8F 

DIMENSION  AUPHA2(1000),Bf TA2C1000),NP3PTS(7,a,18J 
C0MMnN/FIT/IC0EFS,NC0EFS,IEIT(9) 

COMMON/PLPR/MGAM8(7#a) 

RAD»57, 29578 

RE  AD (5, 100)  XSCMIN#XSCMAX, YSCMIN, YSCMAX 

RE  AO (5, 100)  XREF,  YREF,ZREF*DISTG,ObSMIN,ORSMAX 

10  REAOfSp 100)  THETA  f  PH  I 

RE  AO (5, 103)  WGAM, (GAM (I ) , I« i , NGAM) 

READ(B,103)  NPSI,  (PS(I),I»1 ,NPSI) 

103  FO«maT(TS,SX,(5E1O,0)) 

SUMw«380,*(OHSMAX-nBSMlN) 

IF(SltMW,EQ.OtO)  8(Jmw«1,0 

R£AD(5*  lOt )  NFENCE 

IF ( NFENCF ,EQ(0)  GO  TO  1« 

LINES«B1 
ICOEF  8*0 

DO  1?  1*1, NFENCE 
RE AOf S| 1 01 )  NFPTS(I) 

N«NFPT8(I) 

RE  AD (5, 1 00 )  ((FPTS(I,J»K),K*i,i),J«l,N) 

DO  11  J* 1 #  N 
DO  11  k*i,3 

11  XPTSB(,J,K)*FPTS(I,  J,K) 

CALL  CFIT (XPTSB»N, CUE FH» 1 #  S I G  ) 

IF (LINES, GTt30)  *RI TE ( b , 2 1 0 ) 

IF (L INF  8f GT , 30 )  LlNESaO 

210  FORMAT ( 1 H 1 ) 

"RITE  (<N, 211  )  I 

211  FORMiTIlHO, 'FENCF  M|mbFR',I3> 

<nR l  TE  (6 # 2 1 2) 

?1?  FORMAT(1ho,5X, 'BOUNDARY  POINTS  •  BODY  AXES'  /IH0,2(20X, 'X'# 

1  HX,#v*,l!X,'7ft2X)J 
NP*N 

IF(2*(N/?).LT,n)  np*nP-1 
DO  1  SO  J*1#NP,2 
JR  1  *  J*  1 

150  <*ritF(F>,?13)  (M,(XPTSB(M,K)rK3l,3),M*j,JPi) 

213  F ORmaT ( t h  ,2(11X,'(',T?*')',1X,?(F7.2#5X)#F7.2)) 

IF  f2*(w/?),LT,N)  inrITF  (6,?ia)  N,  (XPT9H  (N,x  )  ,ki  1 , 3) 

214  FORiATMH  ,1lX,'(',l2.')##lX,3(F7,2#bX)) 
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>YRITE(ft,21S)  (J.COEFP(J),  J»1 ,«) 

215  FORMAT(1HO,5X, 'COEFFICIENTS  OF  PITTED  CURVE*/!  HO ,  1  OX,  fl  (  •  (  • ,  12 , 

1  ')',Fl0.ft,SX)) 

*RITE(ft,?lfe)  SIG 

21ft  FO«MAT(1HO,5X, 'SIGMA  ERROR  OF  FlTa',Dl2,U) 

LlMES*LINeS+12*N«N/? 

DO  1?  K«i ,U 

12  COFFHFf I,K)sCOEFP(K) 
la  CONTINUE 

15  RE  AO  (5  *  1  0 1  )  NPANEL,lANGLt  IPLOT, IPROB.NDISP 
IFf'JOTSP.EO.l)  M>ISP«<4 
IF(NOISP,EG,0)  NDISPal 

101  FORMAT (1215) 

MAS*0 

00  2000  T  « 1 , N P A N F L 

RE  AO (5. 101 )  I  TYPE,*  PTS,lrOEF3,NBET*,lNTERN 
pEaD(S,  1  00  3  C (XPTS«(l i J) ,  J*1 #3) #L*1 ,NPTS) 

IF  (  I  TYPE  ,EM, 1 )  GO  TO  2fe 
IFfirOtFS.FO.n  GO  TO  20 
IFCir.OFFS.EO, 2)  GO  TO  25 
NCOEFS*! 

00  17  Lai  .ft 
KSL 

IPCL.GT.n  kbL  +  S 
17  IF  I T (L ) sK 

GO  TO  25 

20  PfcAOfS,  1  rip)  (CnFFHCL).Lai.lO) 

MCOEPSsR 

102  FORMAT (ft 0 lo.O) 

25  CONTI WOE 

IF(rcOFFS,fcU,2)  REAOC5.101  )  NCOEFS, (Irl T IL) ,L«1 ,NcOEFS) 

2ft  CA|.l  CFITCXPTSR, OPTS, COfF«,ITYPt, SIG) 

SlGi«t,*SIfi 

IF  (IT  YPE  ,E«,  2,  AND,  *'HF  T  r,T  ,0)  CALL  ADnPTS(XPTSH,NpTS,COEFR, 
l  SIK.MPETw) 

"RITE (ft, 217)  I.ITYPE 

217  FOMMftTClHl ,*WEFLFCTIOM  Panel  NUMBER*, I3,20X, *TYPE*# 12) 

«*RITE  (ft, 212) 

NPeNPTS 

IF ( 2* (NPT  S/2 ) ,L  T , OPTS )  NPsNP«t 

DO  IftO  Jsl,NP,2 

JP1=J+1 

lftfl  ivRITF  (ft, 213)  (M,(XPTSP(M,K),Ksl,i),Msj,jPi) 

IP  (2*(f  PTS/2),LT,NPTS)  tA'RITF  (ft,2l«)  WPTS ,  ( XPTSB  (NPTS ,  K  ) ,  K*  1  ,  1) 

I F ( If OF  E  S , f  0 , 1 )  *  R! TF  f  ft , 2 1  9 ) 

219  FORMaTC Imo.SX, 'COEFFICIENTS  FOR  FITTED  CURVE  ARF  BEING  READ-IN*) 
v"  RT  TF  (ft,21  5)  (J,f  OEFBC  J)  ,  Jsl  ,<i) 

TffITYPE ,E0, 2)  ,'R T Tt f ft , 2 1 8 )  ( J , CUEF H CJ ) , JsS,  \  o ) 

21S  FORMatMh  ,  1  0X,<4f  *  (  *,  12,  *  )  »,F8,R,  1  OX)/lH  ,  1  Oy ,  2  (  *  (  *  ,  I  ?,  *  )  *  ,  FA  ,  «  , 
1  1  0  X  )  ) 

''-RITE  (ft, 21ft)  SIG 

IF  (  i'jF  fc’A'CE.  ,  F  0, 0  ,  OR  ,  I  T  YPF  ,  ftp  ,  2 )  GO  TO  io 


"t? 
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CALL  CHECKS (X,Y,Z,Xl,Yi,Zl,INTC5,l> 

ICOEFSaO 

CALL  CF!T(XPTSP,HPTS,C0EFR0, 1»SIG) 

C  w«ITF(6.«J?)  (CGtF80(LL),LLai#«) 

C«32  FOPMATflHO, 4020,4) 

30  CONTINUE 

00  2000  IQ AMs  1 # NG AM 
^GAMS(I»IGAM)aMAS 
LINF.S«S1 
GAM,<»A«r,Ai«l(  IGAM) 

00  moo  IHSI  =  1,MPSI 
NPSPTSd*  IGA*,  IPS!)«0 
PSIa  PS(IPSI) 

SNPHb8TN(PHI/PAD) 

CSPH*rnS(PHl/RAO) 

S iviTMaSTNl  Tht  TA/WAO) 

C8TH«f;ns(Twf  TA/RAO) 

8NPSaSTN(PSI/PA0) 

CSPSbC08(PSI/RA0) 

DlftCOSf 1 , l)aCSTH*C8HS 
OHCOSn  ,2)«CSTh*SnPS 
niPCOSCl ,3)*-3NTh 

HIpCOS (2i 1 )s8NPH*8mTh*CSP8»SNP8*C8PH 
^lPCnS(2i2)aSNPH*8NTH*SNPS«>CSPH*CSPS 
0I«C0S(2»3)s8mPh*C8Th 

OIPCOSC i» 1 )aC8PH*SMTH*CSPS+SNPS*SNPH 

OIPCf'S(  <,?)■(. SPH*8mTh*SNPS-SNPH*C8P8 

OTPCnS(3#3)arSPH*LSTH 

CALL  Cv»T  (COEFP#OIPCf'S»COEF  I ,  I  TYPE) 

IP  ( UF CF , EG, 0 ,0&#  I  TYPE  , Nfc ,2)  GU  To  AO 
00  iS  t  si ,  10 
35  CPf  FS(I  )=cnf  FI  (L) 

CALL  Cv/OT  (C0tF«Q,l>IWC03,COEFI,  t ) 
ao  co«mt  wf 

on  i  non  miSPai  ,mo)sp 

GAMv/sOrtM^A 
GA  •*1*0,0 

IF(ndisp.Eo.i)  go  t o  as 

TFfiniSP.LF  ,?)  GAMvsGAMV+(«l)**II)ISP/a, 

IFf  IfilSP.GT,?)  GAMLaf-1)**  10  ISP/4, 

4  5  c •T'  T  T  .F 

GVsGAmv  fl«fiA«V 

SlHC  n*-COS(GAW/PAO  )  *CnS  (GAlvL/WAO  ) 
Sli'l(?)sCnS(GAMV/RAD  )  *SIfv(GAML/*AD  ) 

SOMf 3)sSIwfGAMV/PA0  ) 
no  81  LSI  ,0PTS 
00  Si  Is  1  #3 
PP8I ft ,t  ,  J)sO,0 
on  SI  K  s 1  ,3 

51  PPSIU,l,.nsPPSI<l,L,  J)+DIPCOS(K,J)*XPTSR(l,K) 
viRPLS=1 
j  w  P  T  5 ( 1 JsOPTS 
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IF  (NFE>CF  ,EQ,0 )  (,n  TO  5« 

IF  ( I T YPfc ,f  0 , ?)  CALL  THAN(3(RP3I#l#NPTS,COEFI,SUN»n 
00  53  l*t,NFENCE 
00  510  Ja|,« 

510  C0EFD(J)»C0EF8F{L#J) 

CALL  CvHI(CO£FO,OIRCo3,COEFlFri) 

CALL  PI ANPT (COEFIF ( 1 ) #COEF IF (2) » COEFIF (3 ) »COEF IF (8 ) #  SUN ( 1 ) » SUN (2) # 
1  SUN ( 3)  » DISC #0,0I>00 #0,0000# 0,0000 #X,Y,Z1 
IF(0ABS(0ISC),LT,TOL)  go  TO  53 
NFPS*MFPTS(L) 

00  5?  N*1,MFPS 
00  52  J*l,3 
FPSI (N# J)a0,0 
00  52  K«l,3 

52  FPSI (NiJ)aFPSl(N,J)*DI«COS(K,J)*FPTS<t,N,K; 

CALL  SHADOW CBPSI,NPPLS,NPPTS#FPSI,NFPS, SUN, COEF I, GAMV.GAML) 

IF (NWPlS.EO.S)  GO  TO  535 
IF(MPPLS,E«,01  GO  to  poo 

53  CONTI MUF 

535  IFdTVPE.FO.l)  GO  TO  5« 

DO  5a0  TRPL®1 # NPPLS 
LPTS«NBPTS(IRPL) 

CALL  T«AN0(RPSI,IHPL,LPTS#C0EFS#8UN,2) 

NRPTS(IRPL)»LPTS 
500  CONTINUE 

00  505  L«1»10 
5«5  C0FFICL)«CnEF.3(L) 

5«  00  850  IPPLU.NRPLS 

LPT3«NRPTS(IBPU 
DO  800  K»1,LPTS 
MASaMAS+1 
Xi«RP3I (IRPL#K#1) 

YtaRPST(IRPL,K,2) 

Zl*PPSI Cl PPL #K,  3) 

IF ( I T  YPE , EQ , 2 )  GO  TO  56 
00  55  L*1 #  3 

55  VNOBM(U«C0fFI(L) 

GO  TO  58 

56  VNORM( 1  )«2, *COEF I ( 1 ) *X t ♦COEF I  (<I)*Y1  ♦COEFI (5)*Z t ♦COEFI (7) 
VN0PM(2)«2,*C0FFI  (?)*YltCOEFI(«)*Xl*COEFI(6)*ZUCOEF!(B) 
VNORM(3)»2.*COEFI(3)*Zl*COEFI(5)*XltCOEFI{6)*YUCOEF!(9) 

58  ASN»0,0 

ASSNSOaO.O 
00  70  L«l,3 

ASN»ASN4SUN(L)*VN0RM(t) 

70  A83NSQ»ABSNSQ*VN0PM(U*VN0PM(U 

IF(INTEPN,EQ,1,AN0,A3N.GT,0,0)  GO  TO  85 
A3N« A3N/ABSNSG 

IF(AB3(ASN),LT,1,0E»08)  A3N«0,0 
00  80  t*t»3 

90  REFLTN(L)i3UN(L).2,0*ASN*VNOPM(L) 

ICLEAHbO 


39 


IF  (I’.'Ff  *  CF  GO  TO  90 

CAl  (.  H-  TERF  (FPTS,r>!Rr03,,'FFNCF,NFPTS,Xlf  Yl.Zl  ,KFFLTN,ICLFAR) 
iKinjJo.N.ft)  r;n  to  ro 
MS  ALPHAS  1 OOOO, 
nt TAsl onno, 
go  ro  qh 
oo  COMTIWI'F 

tF(DAHS(REFLTN(;i)),GT. 0,9999)  GO  Tn  97 

alphas  t)ATAK.?(-»FFLT  i(5)  #r»S««T  (WF>LTN(  1  )**2+WEF|  TN(?)**2) ) 
hF'TAs  nATAf'2  (9EFLTN(2)  ,  WtFL  T  J  ( 1  )  ) 

9S  alphAskAO  *alpha 
BtTAs  PAO  *BF  T  A 
OG?s-DISTG*OISTG 

call  fin  a  dpt  ( ]  ,000  0,  i ,  noon,  o.onoo,  o.oDnn.n,  onno,  otnooo,o,or>oo, 
x  o,oooo#ntor>oo,nn?#wmTNn ) , 

1  RFFl.Tr  (2)  ,RFFI.  T'  H)  ,Xl  ,  Yl  ,Z1  ,X,  Y,Z.  XI2,  YI2.ZT2,  INTO 
r f f  iMTr,MJ,i )  go  to  <jb 

CALL  3nwM)T(PfFLTNn),PtFL1N(?),PEFl  TN(  *),Xl,Yl,Zj  ,X,Y,Z,  ISHAO) 
IFdSHAO.FN.n  GO  TO  Qh 
X  =  X  I  2 

YSVJ2 

Z*7T2 

Qb  ALPHA2fMA.S)sWAD*liATA'i2("/  +  ZPFF.t)8lJHT(  (  X-XWE  F  )  *  *2*  (  V- YRF  F  )  **2  n 
BETA2('-AS)sOATAM?(Y-YBF>  ,X«XBfcF  )*HAO 
l  F  n  A  0  r;  L  ,  F.  o  ,  1  ,  A  N  0  ,  B  E  T  A  ,  L  T  ,  0 , 0  )  HF.TA*3foO,+HFTA 

if  f  IAiYGL.F'O,  1  ,  ASr'.PFTA2(HAS)  ,LT,  0.0)  MF:  T  A2  f.viAS  )  =  360  .  +  RF  T  A?  (MAS ) 

GO  TO  79S 

Q7  ALPhas-PT  Ft  T1''(3)*90,/DAHS(BFFLTM(  j)  ) 

HE  T  A  =  0 , 0 

PM  ALPHA? (MAS)  ■  ALPHA 
RFTA2('AS)=RFT A 

795  XI  Bsf)IPCOS(l ,  n*XUOI«COS(  1 , 2 ) * V  1  ♦nl«CUS (1 , 3) *Z1 
Ylfl«r>IPC0S(2,  t  )*X)+OT«CUS(?#2)*Y1+OI»COS(2.1)*Z1 
ZIHal)I»C08(?.  t  )*XWDlBf:0sO.2)*Yl*DlKC08n,5)*21 

CALI.  CHf  CKa(XPTS«,NPTS,  X]R,  YiH,ZiB,3I(;3) 

IF (LIhFS.LT ,50)  GO  TO  79*, 

mPTTF(M,200  )  I.GAMMA,  XHfcF, THETA,  Y«F.F, PHI, ZPEF,  OISTG 
?00  F0»mAT(1H1,«0X,«S UW  GLINT  SIGNATURE'  FOR  SURFACE  NO,', 13/ 

X  1H0,SX,»SUM  F.LEVATinM*,Fl?,2,76X,*X«Ff',F9,?/lH  ,5X, 

1  *  A  /  C  PITCH  ATTITUDE', F7,?,76X,'YWEF',F9,2/1M  ,5*. 

2  *  A  /  C.  WOU  ATTITUDE', FM, ?, 76X, *ZHEF',F9,?/1H  ,  106X,  'XOISTG'  ,F7,2) 
*RlTF (6,202) 

202  F0rmaT(IH0|5X.  '  A/C  YA-v',MX, 'SUN  DI SPE RS JON' , 8X , 

l  'ROM  (I  P  A  R  Y  POINTS-SUOY  AXES'*9X»# RE  FLECTION  VECTOR', 16X, 

?  *GL I  NT '/ l H  ,?1X,'VFRT',SX,'LAT',13X,'X',7X,'Y',7X,'Z',I?X, 
i  'ELEVATION'. 3X,  »AZIvUTh', MX, 'ELEVATION ',3X, 'AZIMUTH') 

LINE  9  =  9 

796  nRlTE (6,201)  PSI , Gv , G AML , X 1 H, Y 1 B , Z1 B, ALPHA , BE TA , ALPHA? (MAS) , 

1  BETA2(maS) 

L I NE  S*l.  I WE  8*2 

?0j  FORMAT (1H0,  «X,F7,?,9X,F5,2,3X,F5,?,BX,3F8,?,1 1X,F6.?,4X,F7,?,10X, 
1  F6,?,AX,F7.2) 
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.-a»>'V’»t.«f  ^':^i 


-rm**n>T*  • 


800 

850 

900 

950 

1000 

2000 


2020 

2050 


1  no 


CONTINUE 

NPSPTSd,  IGAN.lPSDBNPSPTSd^uAM,  IPS!)  ♦LPTS 

CONTINUE 

GO  TO  1000 

CONTINUE 

00  950  KS1,2 

MA3«MAS+1 

ALPHA2(mA8)«10000, 

BETA2(NAS)«10000, 

NPSPTSf  I ,  IGAM,  lPSI)«NPSPTSd,  I  GAM,  IPSI)+2 
CONTINUE 

IE(NPANEL,fO,n  GO  TO  2000 

IP(lPLnt,EQ,n  CALL  PL0T{8ETA?,ALPHA2,GAMMA#1#NPSI,NI)I8P, 

1  NP3PTS,XSCNIN,XSCMAK, Y8CNIN, YSCMAX,PHl,THETAtIGAM,I#IANGL) 
CONTINUE 

IF(lPWUHtF.O.O,ANO,IPLOT,PO,0)  GO  TO  2050 

00  2020  I G4N« j ( ngam 

GAMMAsGAM(lGAM) 

IFCIPLOT.EO. lJ  CALL  PLOT (BE T A2 , ALPHA2, G AMMA , 1 , NPSl , NDI8P , NP3PT8 # 
1  XSC*IN,  XSCMAX,YSCNIN#  YS0AX,PHI,TUETA#  IGAM,NiPANEL»  IANGU 
TF(IP«nH,EO,01  GO  TO  2020 

CALL  P«(iRL(ACPHA2,PF:TA2,GANMA,NPSI,NDlSP,NP3PTS,0BSMIN, 

1  OBS^AX.SUnw, IGAn.nPanEL) 

CONTINUE 

PEAf)(5,101)  ICONT 
IFflCOUT.NE.l)  STOP 

GO  TO  10 
FQHNAT (  fe  F 10.0) 

end 
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SUBROUTINE  CHECK«(XPT8B#NPT8,X1B,Y1B#Z1B#3IG3) 

DIMENSION  XPTSS(30#3) 

IFCSIGj.tT, l.OE-OR)  9IG3«0,0001 

DISCS»SIG3 

DO  10  I«1,MPTS 

DI8T«8DRT(XPTSB(I#1)**2*XPT8B{I#2)**2*XPT3B(I.J)**2) 

DIST2«SQRTC  (XPTSB(l,n-XlB)**2*(XpTSB(I,2)«VtB)**2T(XPTSBCl#3)» 
1  Z 1 B 1 **2) 

OISCaOIST2/(Dl3T*91G3) 

IF (0T3C,GT,8IG3I  GO  TO  10 
IF<0ISC,(.fc,Ul8CS)  GO  TO  10 
X 1 M«XPTSB ( I » 1 ) 

V l BaXPTSB ( I #  2 ) 

Z1B«XPT3H(I,3) 

DI3CB«0I8C 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  CFIT(XPTSB#NPTS#COEFB,ITVPE,SIG) 

OIMENSION  XPTSB(30,3) 

REAL*B  SIG# OET » TOL 

PEACES  A(9,9),B(9),XCURF{9,30),F(30),COEFB(10)#COEF8(9> 
C0MM0N/FIT/IC0EFS,NC0EFS#IFIT(9) 

OATA  TOL/t, 00-02/ 

JC*0 

IF(ITVPE.EG,2)  GO  TO  32 
00  20  l«l,NRT8 
F(t)«1.0 
00  20  J*1,J 

20  XCURF(J#L)«XPTSH(L,J) 

CALL  CURFIT(F.XCURF,3»NPT8,C0EFB#A#B#0ET,8IG) 

IFdCOEFS.EO.l)  RETURN 

C0FFR(a)*.0ET 

RETURN 

32  00  35  l=!,NPTS 

XCURF(!,L)«XPT8Ba,  l)**2 

XCURF(2»L)*XPTSB(Lf2)**2 

XCURF(3#L)«XPTSB(L,3)**2 

XCURF(a,L)«XPTSBa»l)*XPTSR(L»2) 

xCiiRF(5,L)»xPTSfl(L,n*XPTSR(L.i) 

XCUMF(R,L)*XPTSB(L#2)*XPTSR(L#3) 

XCuMF(7»l,)«XPT8Ba#l) 

XCURF(A,L)«XPT8B(L#2) 

XCii«F(9,L)«XPTSB(L,3) 

J5  F(L)»1.0 

IFCNC(lEF8,E«,9)  GO  TO  960 

«a  1)0  uh  T*),McnEFS 

A* IF  I T  ( I  ) 

IFCK.EO.I)  GO  TO  96 
00  95  L»1 # n P T 8 

95  XCURF (IfL)«XCURF(X,L) 

96  CONTINUE 

960  CALL  CURFIT(F,XCURF#NC0EFS,NPTS,C0fcFB#A,B,DET,8IG) 
:F(icoFFs,tu,n  return 

IFCNCCIEF8,E0,9)  GO  TO  50 
IF ( IC0FF8#EQ#2)  GO  TO  97 

lF(NPT8tLT,(NCOEF843),OR,3IG,l.T,TOL)  GO  TO  "7 
NCOEFSsNCOFFS+3 
IF(NC0FF8,EQ,6)  GO  TO  99 
GO  TO  32 

97  00  90  I*1,NC0EF3 
9ft  COEFS(I)bCOEFR(X) 

00  90  T«1 ,9 
90  COEFB ( TJsO.O 

00  99  I«1 f MCOEFS 
K*IFITfI) 

C0EFB(k)«C0EF8(I) 

99  CONTINUE 
50  COEFB(10)«-DET 
RETURN 

ENI) 
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SUBROUTINE  CURFIT(F,X,N,m,C0EF,A,B»DET»8IG) 
IMPLICIT  RE AL*8  (A-H.O-Z) 

DIMENSION  X(9,305#Ft30),COFF(10),A(9,9),8(9) 
COMMON/FI T/IC0EFS.MC0EF8# IF  IT (9) 
IFflCOEFS.EO, 1 5  GO  TO  25 
DO  10  I a 1 , N 
8(15*0,0 
DO  10  J* 1 » N 
10  A ( I  ,  J ) *0 , 0 
DO  20  T *  1 » N 
DO  20  J*1,N 
8(J)sR(J)+F  m*X(J,tl 
DO  20  Lsl#N 

20  A(L»J)bA(L»J)*X(J,I)*X(L,I) 

CALL  EDNSOI.  (A,B#N»COEF#DFT,t) 

25  SIG*0,0 

IF  ( ICOFFS.Ea,  1)  nETB-COEF(lO) 

DO  40  1*1 ,m 
SUMSO.O 

no  30  J * l , n 

30  SUMsSMf'+CUFF  f«J)*X(J#I) 

SUM*SUM*DET 
40  SIG*SIGfSUM*SUM 
SIG*DSD«T(SIG/v) 

RETURN 

END 
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SUMHOUT  I  f-'F  KQNSOL  (A,H,I«,X#DET,  NHOMO ) 

IMPLICIT  RFAL*B  (A*H,0*Z) 

OI^ENSTOV  A(9,«n  ,R(9).X(10) 

DATA  TOL/1  .OD-nfl/ 

NRANK*r 

no  in  lei#*' 

10  X(I)«0,0 

t«o 

tCOL=0 
12  HIM 

iCOLsirni-M 
IP  1 aT ♦  1 

1U  AM A  X I »  A ( J , T  COL ) 

L*  l 

IF(IP1,GT#*)  GO  TO  150 
00  15  *aIPl,i'i 

IF (DABS ( A (K, I  COL) ) . G T , DABS ( AmAX I ))  L*k 

15  AMAXiaAfL, ICOL) 

ISO  IF(l>AOS(AMAXI),GT,TOL)  GO  TO  16 
X  (  TCOL )  a  1  #0, 

MRANKSMPANK-1 
IFnCOL.EO.N)  GO  TO  <15 
ICOL=ICOLtl 
GO  Tn  1<1 

16  IFCL.fcG.T)  GO  TO  25 
00  20  Jet,*. 

SvAPaA(I,vl) 

ACT,  JJbA(L.J) 

20  A(L, JJaSwAP 
SwAPaBd) 

P(I)aBa) 

B  ( L  )  a  S  i*  A  P 

?5  X0TVaA(I,IC0U 
DO  30  JaICOL,M 
30  A(I,.J)rA(I,J)/XDIV 

eniaHf  n/xim 

IF(IPl,GTtN)  GO  TO  (1? 

00  (10  KaJPl ,N 
XMtJLsA  (  H  ,  ICOL) 

00  35  J*IC0L*N 

35  A(K, J)aA(K, J)-XMUL*A(I, J) 

90  B(K)afl(K)-XMUL*9(I) 

IF(ICOL,LT,N)  GO  TO  12 

<12  IF(NRANKfe:O,M,OR,NHOMO,EG,0)  X  ( I  COL )  «B  ( I  )  /  A  ( I .  ICOL ) 
(15  IJUMPalCOL*! 

00  50  Tap#  N 
k»m**  I  ♦  1 

IF(X(K),F«#0,0)  GO  TO  47 
I  JIlMPa  IJUMP»  1 
GO  TO  50 

(17  IFINRAMK.EQ.M.OR.NHOMO.ER.O)  X(K)a«(K) 

KP 1 aK ♦ 1 
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00  tih  ,I*K  P  t  >  M 
Ufl  jt  (K)sK(Ki.A(K«IJI]MP, 

SO  CONTIM.If 

het  =  i  .0 

IF(«i«J^K,LTtN)  OFTsO.t) 

wt  r  ijpn 

F  -vJO 


i 
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SUBROUTINE  ADOPTS (XPTS8#NPTS,CQEFH,S I G.^BETW) 

REAI.*B  Cf!EFR(10),CnEFPU0) 

REALMS  SIGfSIG? 

Dl^EMSIOM  XPTSBnO,  J),XPTS(iO,  4),XSTriHE(30, 5) 

REAL*R  A2,R?,C?,02,E?,F2,r.2,H2,K2»L2*  A1  ,B1  ,Cl,Dl  ,A3,B3,C3,DS,E3 
1  F3*G3,H3,K3#Li#SA, SH , SC ,1>ABS , FACT , DI SC , DSQRT 
DATA  TOL/1 ,0F-f!«/ 

A2*COEFB ( 1 ) 

82*COEFB(2) 

C2«C0EFR(3) 

D?*COEFB(«) 

E2*CUEM(5) 

F2*C0EFB(6) 

G2»CCEFH(7) 

*2*COEFH(fl) 

k2*COEFH(P) 

L2*C0EFH(10) 

IFLAG-0 

M*0 

00  100  1*1  #  MPTS 
00  IS  K *  1  *  3 

15  XSTUME(M,K)*XPTSH(I  ,K) 

IP1*I+1 

!F(I,E0PWPT8)  IP1*1 
IFflFLAG.EQ, 1)  GO  TO  200 
5  x*(xpTsea.n*xpT8P(iPi#n)/2, 

Y*(XPTSB(I,2)tXPTSB(IPl,?))/2, 

z*(xpTsen,3)*xPTSpnpi#3))/2, 

OfeL*A2#x**2*B2*Y**2+C?*Z**2^02»X* Y^E2*X*7*F2*Y*Z*G2*X+ 

1  H2*Y+K2*Z+L2 

IF(A8SCQtL),LE,3t*SIGl  GO  TO  100 
1P2*TP1+1 

IF(IPl,EQtr!PT8)  IP2*1 
00  20  *«i  #3 
XPr8(l#H)*XPT80(I*K) 

XPTS(2,K)*XPT3R(IPl ,K) 

20  XPTS(3,K)«XPTS8(IP2*x) 

CALL  CF!T(XPT3,S,C0EFP,1,SIG2) 
f l*COEFP ( 1 ) 

Bl*C0EFP(2) 

C1*C0EFP(3) 

01*C0EFP(«) 

GO  TO  208 

200  OEL*Ai*XPTSB(IPl , 1 ) ♦» 1 *XPT30 ( IP t » 2) *C 1 *XPTSB ( IP  1 , 3) *01 
IF(A88(0EU.LE,ToL)  GO  TO  20« 

IFLAGbO 
GO  TO  5 

20«  00  205  K«1,J 

XPT8(1,K)«XPTS8(I#X) 

205  XPT8(2,K)*XPT8R(IPl ,K) 

208  IF(DAB8(Bn,GT,0A88(Al),AN0,DABS(Bl).GE,DA8S(Cn)  GO  TO  30 
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if  (iiAMsrc  i ) ,(-r  ,oflHS(  ai  )  ,flf-n,t.)AHS(cn  .GT.DAosfBi ) )  fin  to  «o 
IKIFlAfi.fcn.n  fin  TO  pop 

r.3*C2*«?*(n/An**?-f  p*ci/ai 
f  3*F ?♦/,* A?*r  i  *hi  /  a  i  **p-np*ci /ai-fp*«i /a  i 
H5sH?  +  ,?,*A?*ni*Hi/Ai  **p-n?*ni/Ai-G2*Rj/Ai 
■<3aK?  +  ?,  *A?*i'l  *C1  /AlW*F.?*01/Al-G2*Cl/A1 
l3slp*ap*(i>i  /a j  'M/ai 

?09  IF  (  APS(XPTS(P,2)..XPTS(  1  ,?))  ,GT.  ARS(XPTS(2,3}*XPTSC1  ,  3))  )  SO  TO  2« 

no  P30  js l ,  luHF Tiai 

VSM+1 

ZsxPTSd  #3)>PLOAT(^)*(XPTSf^,3)«XPTS(|,3))/PLOAT(AiRFTiW*l  ) 

SAsrt3 

S*sF 3*?+h3 

SCsC3*?*7^«3*Z+L3 

IF  (OAH.S(SA)  ,LT,  1  ,on-OH)  fin  to  2  3 

FAfTsl.O 

0ISCs8H**?*«t *SA*SC 

IFfOTSr.l  T.o.o)  «HI  Tfc  IB. 600)  OISC 

IF (DISC .IT, 0.0)  DlSCsO.O 

*,00  FOPwaT  (///,  ts(  '**)  ,5*,  #n  1  SC«  I  ^  INAMT  is  lfss  than  zfpp  in  adding  ro 
IiivoAkY  POINTS*  ZF  Wf1  ASSUfcK>#/lHO,?OX,  'ACTUAL  VALUE  T  $  •  ,  D 1  A ,  a) 

Pjo  ys-Srt/^,/SA  +  FACT*l)S"PT(r  !SC)/2,/SA 

DISTsARSfXPTSn  ,?)-Y)+AHS(XPTS{2#2)*Y) 

IF  <  r  ,1  T#XPT8n#2)#AN|.,V,LTtXPTS(i,2J)  GO  TO  2? 

IF  f  Y.fiT.  XPTSd  ,P),A.'ir.,Y,GT,XPTS(8»?))  GO  TO  ?.? 

P15  X  STUFF (M|?)SY 

XSTOWf.  f  '^Oa-Dl/Al-Cl  *Z/A1.Ht*Y/Al 
XSTOhE ( F  ,  3) »Z 
GO  rn  PSO 

??.  IF  (F  ACT  ,  L  T ,  ft  ,  0  )  GO  TO  220 
FACTs-1  ,0 
OIST1  af'IST 
GO  TO  ?10 

PPO  IF (DTST.LT.DlSTi)  GO  To  PtS 

YS-SH/2, /SAfDSOPT (DISC )/?,/SA 
fin  TO  PIS 

P3  IF (OAHS(SB).l T.TOL)  GO  TO  500 
ys-SC/S^ 

GO  To  p 1 5 
230  CONTINUE 
GO  TO  00 

pa  00  2B  a  *)  «NHF.  Ti'1 
msm+1 

YSXPTS n,2)*FLOAT{N)*(xPTS(?,2)-XPTS( 1,2)} /FLOAT (NBfTw^l) 

S  ArC  3 

SO  =  F  3*Y*K3 
SC«B3*y*Yfh3*y*L3 
IF{0ABS(SA),LT,1 ,00-ftR)  GO  TO  IS 
FACT«1, 

0ISC«8P**2»fl,*8A*SC 

IF(f)ISC,t  T.O.rt)  "HHElAiFOnj  DISC 
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»«^-»  %r  wv**.-v  ,*..  -*  ^ar.  .,,•  v.  *«*<«,  .-,  .-  ...„  vV..e,«~  .  ^  ^..v.vr*r  *r,Y'*  .«-*-«■*•  .  v^w^T  "*:**' •*-'• ^ 


IF (DISC.LT ,0,0)  DISCiO.O 
240  Z«-8B/2,/3A+FAiCT*D3QRTCDI8C)/2,/3A 

8l3jr»AB3CXPT3<l#*l*Z)*AB9CXPT8C2»S)-Z) 

IF (Z,LT,XPT3( 1 , 3) , AN0,Z,LT,XPTS(2, 3) )  GO  TO  242 
IFCZ,GT,XPTS(1#3),AND,Z.GT,XPT3C2,3))  GO  TO  242 
-m  XSTOREfWTH*Z 
X8T0RE(M,2)hV 

X9T0RE(M,l)«-01/Al.Cl*Z/Al.Bl*Y/Ai 
GO  TO  20 

242  IF (FACT ,LT ,0,0)  GO  TO  245 
FACTa-1. 

DiSTt*Ol8T 
GO  TO  240 

245  IF(DI9T,LT,DI8Tl)  GO  TO  241 

Z«*8B/P,/8A*030RT(0!9C)/2,/8A 
GO  TO  241 

25  IFCDAB3CSB),LT,TOL)  GO  TO  500 
Za.SC/8B 
GO  TO  241 
28  CONTINUE 
GO  TO  90 

30  IFCIFL4G.E8.il  GO  TO  309 

A3aA2*B2*(Al/Bl)**2-02*Al/Bi 

C3aC2+82*(Cl/fll)»*2«F2*Cl/Bl 

E3aE2*2t*H2*Al*Cl/BU*2«02*Cl/Bt«F2*Al/Bl 

G3aG2+P,*B2*0l*Al/81**2-O2*Ol/Bi«H2*Al/Bl 
K3aK2+2,*B2*Dl*Cl  /Bl**2*F2*Dl/Bl.H2*Cl/Bl 
L3aL24«2*(0l/8n**2-H2*0i/Bl 

309  IF(A0S(XPTS(2,l).XPTS(t,l ) > ,GT , ABS C XPTS C2 , 3 J-XPTS C 1 # 3) ) )  GO  TO  34 
00  330  N*1 , NB6TW 

M*M+1 

ZaXPTS ( 1  *  3 ) 9FL0AT (N )* ( XPT8 (2»  3)»XPT8( 1 . 3) ) /FLOAT (NBETW^l ) 

3AaA3 

3Bafc'3*Z+G3 
3CaC3*Z*Z+K3*Z+L3 
IF(OABS(SA),LT, 1,00-00)  GO  TO  33 
FACTal, 

0ISCaSB*SB-4,*SA*SC 
IFCniSC.LT.O.O)  WRITF(0»6OO)  DISC 
IFC0I3C.it. 0,0)  01  SC *0,0 

310  Xs-3B/?,/SA*FACT*DSQRT(DISC)/2,/SA 
0I3T*AB8(XPT8U,  .  )-X)+A88(XPT8(2,1)-X) 

IFCX.LT, XPTSCl.l), AND, X,LT,XPT3(2#1))  GO  TO  32 
IFCX,6T,XPT8C1#1).AND.X.6T,XPT8(2»I>>  GO  TO  32 

315  X3T0RE(M,l)aX 

X3 TORE  CM, 2) *-01  /B 1 -A  1 *X/B1 -C 1 *Z/B1 
XSTORE (M, 3)«Z 
GO  TO  330 

32  IFCFACT.LT, 0.0)  GO  TO  320 
FACTa-t, 

OISTlaOIST 
GO  TO  310 


i 

1 
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£*Jw5*WSM(a*Sli82? 


320  IF COIST.LT.OXSTi )  CO  TO  315 

X*.SB/2./8A*03QRTCOISC)/2,/3A 
GO  TO  315 

33  IF(OABS(SB),LT,TOL)  GO  TO  500 
X*-SC/SH 
GO  TO  315 
330  CONTINUE 
GO  TO  90 

3«  00  38  M«1,NBETW 

MsM+] 

X*XPTS(1,1)*FL0AT(N)*(XPT8(2,1)*XPTS(1»1))/FL0AT(NBETW*1) 

3A«C3 

SB«E3*X*K3 
8C*A3*X*X+G3*X*L3 
IF(0ABS(3A),LT. 1,00-08)  GO  TO  33 
FACT* 1  , 

DISC* 3 8*88-4, *SA* SC 

IFC0I3C.LT, 0,0)  '"RITE  (6.600)  DISC 

IFC0I3C.LT, 0,0)  nisc*o,o 

340  Z«-88/2,/SA+FACT*080RTC0ISC)/2,/3A 
0IST*A88(XPT3Cli3)-Z)+*PS(XPTSC2»3)-Z) 
IF(Z,LT,XPTS(l,3),AN0.Z,LT,XPT3C2,3))  GO  TO  342 
IFCZ,GT,XPTSC1,3),ANP,Z.GT,XPTS'C2,3))  GO  TO  342 

341  XST0RE(M,1)*X 
XST0RE(M,2)*-ni/Rl«Al*X/81»Cl*Z/Bl 
XSTORF (M»3)*Z 

GO  TO  38 

342  IFCFACT.LT.O.O)  GO  TO  3«3 
FACTb-1  , 

0I3T1*0I8T 
GO  TO  340 

343  IF(0I3T,LT,0!8Tl)  GO  TO  34J 
Z*-SH/2,/SA*03f3RTC0ISC)/2,/3A 
GO  TO  341 

35  IF(DA8SCSR),LT,T0L)  GO  TO  500 
Zs-SC/SB 
GO  TO  341 
id  CONTINUE 
GO  TO  90 

40  IF(IFLAG,EO,l)  GO  TO  40R 

MxA?+C2*(Al/Cl  )**?-E2*Al/Cl 

83*R2+t2*(Bl/Cl )**2-F2*Bl/Cl 

03*02*2, *C2* A l*Bl/Cl**2-F2*ftl/Ci-F2*At /Cl 

G3*G2*?,*C2*D1*A1/C1**2*E2*D1/C1-K2*AI/C1 

H3»H2*2,*C2*01*Bl/Cl**2"F2*0l/Cl-K2*Bl/Cl 

L3*L2*C2*CI>1/C1  )**2»*2*U1/Cl 

409  IF C AHSC XPT 8 C2»?)-XPTSC 1 »2) ) ,GT , ABSCXPT3C2. 1 )-XPT5C 1 » 1 ) ) )  GO  TO  44 
00  430  N*1 ,N0f  TW 
M***1 

XBXPTSC1, 1) ♦FLOAT (N)*CXPTSC2,1 )-XPTSC 1 , 1 )) /FLOAT (NBE TW* 1 ) 

SA*83 

SB*03*X*H3 
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4»H  **n»  <w»-  K-x»,r»»"»<*«e»t«-’  •  tt  »<"•»<«»**  r-'JwwVf- * 


»■ 


3C«A3*X*X*G3*X4L3 
IF(DAB8C8A),tT,l,00o08)  GO  TO  45 
FACT«1. 

0ISC*SR*SB.4,*SA*SC 
IFCOISC.LT, 0.0)  WRITEC6,600)  DISC 
IFCOISC.lt, 0.0)  DISCsO.O 
ato  Y»*SR/?,/SA*FACT*0SQ«T(DISC)/2,/SA 

DIST*AflS(XPTS(l ,  2 ) «  Y )  «.ABS  C  XPT8  (2 , 2 )-Y ) 
IFCY,LT,XPTSn,2).AND.Y,LT.XPTS(2#a))  GO  TO  42 
IF CY,GT.XPTS(1»2).AN0,Y,GT,XPTS(2#2))  GO  TO  42 

aiS  XSTOK£(M, i )sX 

XSTURt:  (M#2)»Y 

X3TORF:c<VS)«-01/C1-A1*x/C1»81*Y/C1 
GO  TO  43« 

42  IFCFACT.lt, 0.0)  GO  TO  420 
FACTs-1  , 

0IST1«0IST 
60  TO  410 

420  IF  (OIST.LT.OISTl)  GO  TO  415 

Ys.SK/2./SAtnSOR7 (DISC) /2,/SA 
Gn  TO  415 

45  IF COAHSCSB) .LT.TPL)  60  TO  500 
Ys-SC/SB 
GU  TO  415 
430  CONTINUE 
GO  TO  40 

44  00  4*  ^■l.MBFTfc 

YaXPTSM ,2)fFL0ATCF )*(XPTSC2,?)*XPTSCl ,2) J/FLOAT  CMBET^^l ) 
S  A  »  A  3 

SHa03*Y+G3 

SC*B3*Y*YtH3*YtL  3 

IF COARS(SA) ,LT, 1 .00*08)  60  TO  4b 

FACT*! , 

ni$C*SR«SR*4,*8A*SC 

IFCOISC.LT, 0,0)  ^«TTEf6»600)  DISC 

IF  (DISC, IT, 0,0)  DISJCsO.O 

440  Xs*SR/?,/SA+FACT*DSO«T(UISC)/2,/SA 
OISTsAPSCXPTSCI . 1 )-X)+AflS(XPTSC*# 1 )*X) 

IFCX.LT, *PTSn»l).ANi\X#l  T,XPTSC2,1))  GO  TO  442 
IF(X.GT,xPTSCl,l),AM7f X , GT , xPTS (2 , 1) )  GO  TO  442 

441  XS70RE  (5, j ) «X 
XSTORE(V,2)*Y 

XSTORE (^,3)**D1/C1*A1 *X/C1*HI *Y/C1 
GO  TO  48 

44 2  IFCFACT.lt, 0,0)  GU  TO  443 
F  A  C  T  e  •  1 , 

01  ST  iso  I  ST 
GO  TO  440 

443  IFCOIST ,L  T.OISTl )  GO  TO  44J 
Xs-S8/2,/SA^DS«HT(0ISC)/2,/SA 
60  TO  W41 
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«5  IF(0AB3(S8) ,LT,TOL)  GO  TO  500 

X«-3C/SB 
GO  TCi  441 
48  CONTINUE 
RO  IFLAG*1 
100  CONTINUE 
NPTS«M 

00  HO  I»1#NPTS 

00  110  jst,5 

110  XPTSH(I,J)aXSTORE(I,J) 

RETURN 

500  *RITEC6,610)  SA,SB,SC 

STOP 

610  F0R'tAT(///,l«>C'*'}>5X#'QUADRATIC  EQUATION  SOLUTION  FAILS  IN  ADDING 

1  BOUNDARY  POINTS-  CHECH  INPUT  •  PROGRAM  WILL  STOP  HERET//1HO,20X, 

?  'VALUES  OF  COEFFICIENTS  OF  QUADRATIC  EQUATION  ARE  A«' , D15 . 4 » 2X , 

2  'Ba'.D!5.4,2X, 'C*'#015,4) 

END 
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SUBROUTINE  CVfll {COFFR,OlKCOS»COfcFI, ITVPF) 

IMPLICIT  REAl*8  (A-H,0«Z) 
n I  WENS  ION  OIRCOS(3,3),COFFI(lO)#CnEFB(10) 

IF(TTVPF.FO.?)  GO  TO  S02 
00  50  I  a  1 , 3 
COFFHUoO.O 
00  SO  jsl ,3 

50  COFFT  (UsCOFFI  (L)  ♦OIRCOSf  J#L)*COfcFH(J) 

C0FF1 CaJsCOfFBU) 

return 

502  COmT  I  mi'E 

00  S06  1*1 .3 
SAaOIRCOSCl ,L) 

S8aOIRCOS(?#U 
SC«0tRcnsc3,L ) 

COf  F I  (L  ) «C0{  F  B(1 )  #SA*8AfC0f  F8 (2)  *8B*8fl+C0tFR ( 3) *SC*Sf-  +  COEF 8 ( <0  * 

1  SA#SH*CnfcFH(S)*SA*Sr+COEFRt(b)*3«*3C 
IF(L,R0.3)  GO  TO  50« 

3A2*r)IRC08(l  * L  + 1  ) 

S0?BOIRCOS(2iL  +  n 
SC2*niSCOS(3,L*n 
GO  TO  505 

50a  SA2a0lHC0S(l ,2) 

S82Bni«C08(2.2) 

SC2*t>lRCns(3,2) 

505  CUEFI  (l  +  3)«2,*COEFB(  n*SA*SA2^2,*C('tFB(2)*SB*Stt2f2,*C0FFB(3)*8C* 

1  SC2  +  rnFFHm*($A*Sn?  +  SA2*SB)+C0FFB(5)*(3A*SC2*9A2*SC)+C0EFB(M* 
1  ( SB*SC2+SB2*Sr ) 

COEF  1(1 +6)80,0 
00  506  K»l,3 

506  COEFia>6)«enF.F!(l*6)+OIRCOS(K,L)*COEFBCK*6> 

COEFI(10)«COEFB(10) 

RETURN 

ENO 
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SUBROUTINE  TPANQ(RPSI,IPPL,NPTS*COEFI,8UN,IC> 

DIMENSION  COEPI(10),8IJN(5)  ,RP3I(5,50,3) 

RE AL*8  COEFI,8UN,A,B,C,D,£,F,G,M,CK,CL,SL,SM,SN,DISC,Xl*Yl,Zl , 
1  X, Y,Z,XQ2, YG2,ZQ2,DSQRT 
DATA  ToL/l,0E-0«/ 

COMMON/SAVE/SL#SM,SN#OISC.CSTH 

AbCOEFI(I) 

B*COFFI (?) 

C*COFF  T ( 3 ) 

0«CO£FT(a) 

GO  TO  (10,505,  IC 
10  SLbSUN(I) 

SMb3UN(2) 

SN«SUN(3) 

C  wBTTF (6,200  5  A , B, C , D , SL , SM, SN 

CSTHs(A*SL+B*SM*C*SN)/OSQPT(A*A+B*8*C*C) 

IF(AHS(CSTH) ,LT,0, 1 7)  RETURN 
DO  20  IH.NPTS 
X 1 sPPSI ( IPPL  *1,15 
Y 1 sRPSI ( IPPL  *1*2) 

ZIbRPSKIPPL,  1,3) 

CALL  PI  ANPT(A,B,C*D,vSLi8m,SN,DISC,X1  ,Y1  ,Zl,X,Y,Z) 

C  aRITF(*,200)  X1,Y1,Z1 ,X,Y,Z,DISC 

C200  FUHMAT(1mo, 7016.il) 

IP  (OARS  (DISC.)  ,LT,  TOL  )  RETURN 
PPSI(1PPL»I»1)BX 
HPSI(IWPL,I#2)BY 
20  RPSI ( IPPL, I,3)»Z 
RETURN 

50  IF(DAHS(DISC) .LT.TOL)  RETURN 
TF(AHS(CSTR).LT,0,17)  RETURN 
EaCOEFI(S) 

FsC0FFI(6) 

GBCOEFK7) 
hbCOEFI (B) 

CKbCOETI  (9) 

CLaCOEFI(lO) 

c  lip  I TE  (6,200 )  A,W,C,D,F,F,G,H,CKfCL,SL#SM,SN 

Lso 

00  50  1*1 ,NPT8 
X 1 sPPSt ( I  PPL ,1,1) 

YlaRRSi (IRPL#  1*2) 

Z 1 bRPSI (IPPL » 1 ,3) 

CALL  OHADPT(A,H,C»n#F,F,G.H,CK,CL*SL,SM,SN,Xl , Y1,ZI#X# V»Z» 

1  XQ?,YO?*ZU?» INTO 
IF(INTr,EQtO)  Y'RITE  (6,500) 

IF  (IMC.FO.O)  GO  TO  50 
P  wRI TE (6,200)  X,Y,Z, XN2* Y«2,ZU2 

IF(INTC.EO.I)  GO  TO  u 5 
CALL  CHECK5(X, Y,Z,x 02, YU2, ZM2, IN TC5* 2) 

IF ( InTCS.EM, 0 )  wu I TE (6 , 500  ) 

IFdMTfS.EO.U)  GO  TO  50 
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IF(INTC5.E8,1)  GO  TO  «5 
L»LM 

«P3I (IRPL,t,l)*X02 
RPSI (IRPl,L,?)sYQ? 
pp3i(ipPL,i#:n*zo2 
P5  L*Ltl 

PP8KIPPI  ,L,n«X 
RP3ICIRPL#L,2)«Y 
RPSI  (I»PL#L,3)«Z 
SO  CONTINUE 
MPTSsL 

500  FO«MAT(1HO,10('*')*»OUE  TP  CURvATURF  OF  SURFACE  *, 

1  'CERTftlN  BOUNDARY  POINTS  fclLL  BE  MISSING  IN  OUTPUT'/ 

?  1H  ,10X,'F0R  THIS  COMBINATION  UF  SUN  ANGLE  AND  A/C  AZIMUTH* ) 
C 

RETURN 

END 
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SUBROUTINE  CHECK5(X,Y,Z,XG)2, YQ2, ZQ2 , I NTC5, IC ) 
COMMON/CTRQ2/OIRCOS(3,3),XPT8B(30,3) ,XMIN(3),XMAX(3),NPTS 
DIMENSION  X8  (  3)  #  XQ2B  (3) 

REALMS  X,Y,Z,X02, Y02,Z«2,DIRC0S 
GO  TO  (10,30),  IC 
10  00  20  1*1,3 

XMAXm«XPT3B(l,I) 

XMTN(I)«XMAX(I) 

00  20  J«2,NPTS 

IP  (XPTSBfJ,  I)  .LT.XMIMI) )  XMjNd  )sXPTSH(J,I ) 

20  IFfXPT8«(J,I),GT,XMAX(I))  XMAX ( I ) bXPTSB ( J , I ) 

RETURN 

30  XB(n«PlRCOS(l,l)*X*OIRCOSn,2)*Y*OIRCOS(l,3)*Z 
XB(a)*0lRC0S<2, 1 )*X*0IRCoS(2#2)*Y*DIRC0S(2,3)*Z 
XB(3)*0I«C0S(3, 1 )*x  +  OI«COS(3,2)*Y*DIRCOS(3,3)*Z 
XQ2B(1  )»f)IRCOS(l  ,1)*x«240IRC0S(  1  » ? ) *YQ2*DIRC0S ( 1 , 3 ) *Z02 
XQ2R(2)al.)TRCoS(2,n*X«2fOTKC0S(2»2)*Y«2+0lRC0S(2»3)*ZO2 
*Q2B(  3)«0IRCPS(  3,  l)*XO2^OIwC()S(3#2)*Yij24>0lRC0S(5#3)*ZR2 
IX*() 

I  X?aO 
INTC5*n 
00  <40  K a  1  ,5 

IF(XR(* )  ,IT# (XMIN(X )-0,0%*(XMAX(K).XMIN(K)  )-0,01))  IXsl 
IF(XB(K)  ,GT.  (XMAx(K)>0,US*fXMAx(X)«XMlN(K)  >4-0.01))  I  X»  1 
IF ( XQ2R ( *  ) ,LT,  ( X*lN (K) -0,05* (XMAX (K)-XMIN(K) J-O.Ol  ) )  1X2*1 
<40  IF  (XG2K(K)  ,GTf  (X^AX(K)+0,05*(XMAX(K)-XMIN(K))+0,01  )  )  1X2*1 
TF(lX,FQ,n  RO  TP  60 
IF  (Ix2.ro, 1)  GO  TO  50 
I  pvl  T  C  5  *  2 
RE  TURN 
SO  iNTCSal 
RETIIRm 

60  TF(IX2,F«,1)  RETURN 

IMTC5S1 
X3X02 
Y*y02 
Z*Z02 
RE  TURN 
Cviii 
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SII'WU  T  I  ME.  PL  AK-PT  (  A  ,  b » c »  n  #  SL  •  8»*»  SM  #  0 1  SC  #  x  1  ,  V  I  ,  Z I  #  y  ,  V ,  7  ) 

IMPLICIT  RfAl,*«  (A-H,(W) 

A8S(XX )sn>APS(XX} 

0  A  r  A  TOL/1  ,0f-04/ 

IF (APS (SI  )  ,r,f  ,  AH5(S^)  ,A.H»,  AHS(!jL)  ,(;e  ,A»S(8^)  ) 

1  C*3''/SL 

(F  (A«S(  St  )  ,LTf  AHS(SM  ,  AMT,  ABS(SM)  ,  L  E  ,  ABS(SN)  ) 
i 

If (APS  (SI  )  .IT,  ABStSM .AMD. AHS(S^) ,LT, A«8(S*)  ) 

1  H  ♦  C  a  S  *  /  S  M 

IF  (AttS(OTSC)  ,LT,TO|.  )  WF  1 « 

IF  (APS  (SI  )  ,LT,AHS(SN)  ,  AMD,  AAS(SM)  ,U' .  ABS(SN)  ) 

IF  (AHSfSl.  )  ,LT,AHS(SH)  ,  AwH,  A8S(S;'»)  .L  T.AHS(SM)  ) 

XS  (-0-P*  Y  1  +  H*X  1  *Sf1/SI-  *c  *z  1  +  c  *X  1  *Sf'l/SL  )  /  DISC 
YSYU8^(X«X1  l/Sl 
ZsZUSf  *(X-XI  )/S|_ 

METI.IRM 

1  0  la (-n-A*xl+A*Zl*SL/SN-P*vt  +  8A/1 ASM/SNJ/OI8C 
Xsxl ♦SI  *(Z-/1 )/SM 
YBY1*8^*(Z-Z1 )/SN 
Rg.  7  iJP‘m 

?n  ys(.|i-A*x  )  ♦  «*SL*y1  /  S^»c*71  ♦CaSMayI  /S'*)  /(A*SL/SM*B>CaSM/S*1 
xsxttsi.*(Y-vn/sf 
Z«ZltSN*(Y»Yi)/S*‘ 

RF  T'JWfJ 
f  MO 


OISCbA+BaSM/SL* 

DlSrsA*SL/SN+ 

0ISC«A*8L/SM> 

bn  to  io 

bO  TO  ?0 


3 

3 


i 
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SU*«0UTI*R  SHADOW (RPai.WRPLS,NRPT3,PP8l»NrPTS, SUN, COEFI,GAMV,GAML) 
01  PENSION  WPS  1(5, 30, 3) ,N«PTS(5)#FPSI ( 30 , 3 ) , SUN ( 3) , 

1  CnEPI(tO),PP3IN(5, 30,3)  ,«PS(30,2),N»PTSN(  5) ,8P8(30,2) , 

?  8PS! (30,3) 

DHENSIOt'  WPSN  (  30  #  P ) 

WfeAt*8  SUN,CO£FI,A,8,C,f>,SL,SM,SN,OI8C,X«,  YO,Z«,XX,YY,ZZ,DSPRT 
DATA  TPL/1.0f-0«/ 

AsCOEFI ( 1 ) 

P»CDFPI(2) 

C=COFFT(3) 

OsCOFFI  (<l) 

SL«SIJN  ( 1 ) 

SMsSlJM  (2  ) 

S'M*Slin(3) 

C3Th*( a*SL+R*Sm*C*SN) /USGWT ( A*A*B*fUC*C  ) 

IKARS(CSTH),LT,0,J7)  RETURN 
N'SsO 

DO  10  J«1,NFPT3 
xa*FPS i (i,i) 

YUaPPSI (1,2) 

ZA«FPST(!,3) 

CALL  PL ANPT(AfK,C,n,SL,SM,SM,0I3C,xa,Ya,Z«.XX, YY,ZZ) 

IF  (DAH.s  (DISC)  .LT.TOL)  KETUWN 

CALL  SDWMJT(SL,3M,SN, X«, YA,Za,XX, YY,ZZ,ISHAD) 

I F  ( I  SH  A()  ,£0,0)  GO  TO  10 
N3aM8+1 

CALL  PL ANPT(SL, SM,SN, 0,0000, SL,SM, S*,OI  SC, X«, Y«,Z«, XX, YY,ZZ) 

SPSI (MS, 1 )SXX 
3  P  3 1 (NS , 2 ) *Y Y 
SPSI (NS* 3)«ZZ 
10  CONTINUE 

IF(MS,LE,2)  RETIJWN 
nwpL8N«0 

00  120  IWPLal  .tMWPLS 
!Ca2 

IF(IHPL,EQ,1 )  IC*  I 
NWiNWPTS ( IWPL ) 

00  15  I«1  ,NR 
X<l«RP8I(IRPL,I,l) 

Y«iWP3rn«PL,  1,2) 

ZR«RP3!(lRPL,T,'l) 

CALL  PLANRT(SL,SM,SN,O,0D00,SL,SM,8N,DJSC,X«,Y4,Z«,XX,YY,ZZ) 

RPSI(IHPL,I,n»XX 

RP3I ( I PPL ,I,2)*YY 

15  BPS  I C I PPL , 1 , 3)*ZZ 

CALL  TPAA/SP(RPSI,IPPL,'VR,SPSI,N3,RPS,8P3,IC,GAMV,GAML) 

IR*1 

16  XC«RPS ( IW ,  1 ) 

YC«RPS(IP,2) 

CALL  CIRCLE(SP8,N8,XC,YC,INCR) 

IF(INCR,NE,2)  GO  TO  17 
tF(XR.EQtNR)  GO  TO  120 
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I  Pa  I R ♦ 1 
GO  TO  lf> 

17  f  F  L  A  r,  s  *) 

IFCtMCB.EiJ.t)  TFIARsr 
Mf  X tso 

IFfTFLAS,  EO.M  MBMaO 
lF(IFLAK,fcU,#.)  BO  TO  25 
MRPLSGxORPI.3M+l 
UR’.la  1 

»PSMf  1  .  Dave 
PPSOfl  ,?)aVC 
I F  t AGs  1 
?5  1 9P 1  a  1 9 ♦ 1 

IF(I«#Mi,Nfc)  T  RP 1  s  i 
IonuPLso 

so  XtsPP3(tP,n 
YlaPPSC IW,P) 

S5  X2aMPS  f I RP  t  #  1  ) 

Y?sKPS(IKP1 ,2) 

CALL  If'T£Pr.  (SPS.MS,Xi  ,>1  ,  X*,  Y2,X,  Y,KL,  IF  AIJLT,  IDDURL, 
1  X4I.SH,  YAl.SO) 

IP  f  IPLAG.fctJ.f),  AOO.IFAUlT.FG.O)  GO  TO  70 
l  F  ( I P  l.  A  (1 1 F  G .  1  t  A  NO ,  I F  A  uL.  T ,  E  0 , 0  )  GO  TO  «0 
I F  ( I F  L.  A  G ,  F  0 1 9  ,  A  NO  #  I F  A  t  iL  T  ,  P  i,)  #  0 )  Go  TO  00 
IF  f  IPLAG.FO.ft)  GO  TO  i8 
CALL  Cf-fcCKl  (Xi»,  Y?#RPS^»NRKif  ICHK) 

IF(TCHK,F«#n  nn  TO  on 
I F  L  A  G  a  1 
N  PMs  NR*" ft 

RP8P»CMPN»1  )  =  *2 
RPS«( NRN,2)»V? 

SB  IRsIRPl 

IF  f  JRPl  ,F(3,1  ,AND,lFl.AGtt(j,6)  GO  TO  90 
TRP1 sTR+1 

T F  f  I W  ,  F (J ,  M R )  I RP  }  a  1 
IOOURLcO 
GO  TO  10 
40  IFLAGs? 

MFXTsIPPI 

XSFPsX 

YSEPeY 

CALL  CHFCK1 (X, Y,RP3N,N«N, ICHK) 

IF(lCH«,£0.n  GO  TO  90 

MR.\j*NRfj*l 
RPSNfNRN, l)aX 
RP3N(NRN,2)aY 
ISaKL 
ISPlalS+l 

IF(IS.fcO.NS)  ISPl«J 

IOnUBLal 

XAlSOax 

YALSOaY 
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KO  f *BO 

so  xisspsns.n 

YlaSPSf IS, 2) 

*P*SPSf ISP1 , 1 ) 

Y^aSPSf  TSP1 ,2T 

CALL  U  T^^C(«P3.NH,m  ,Yl,»?,r?,X,  V.KL,  tF  AUL  T  ,  I UfhlBL  #  X ALSO,  V  AL SO ) 
IP  (IFLAG.fU,?, AMO.TF A(tlT,I0iO)  GO  TO  75 
IP'IIFLAG.eO,  l.ANO.TF  AilLT.EO.O)  GO  TO  75 
IF  UH  AG.F'J,  5,  ANOttfOT»,E0,  l  )  r,0  Til  60 
IP  UFIAG,f,),3)  GG  TO  s* 

Tr(TPLAG,PO,7,»NntIPAiiLT,TG,0)  GU  TO  7 S 
IF(  lH.Afa,Eu,7,  ANI»,M>T«,F.Q,I )  GO  TO  5« 

CALL  ( l«eLF(FP8ffcO,Xt .71,I*CO) 

IPn*'C*-,PO,l  )  GO  TO  S* 

If ( IFLAG.ru.  1)  r,n  to  jo 
54  CALL  CTOCLl  (PP8, *»,<*,  V2.PCH) 

IK  fifw.PQ,  11  GO  TO  60 
IP ( IFLAG.P 0,7 )  I f LAG* 1 
IF (IPLAG.EG.J)  GO  TO  50 
IP  (IFLAG.EG,7)  I  FI.  AG*  t 
t  f  L  A  G  a  7 
I  OOUHt.  *  t 
YAI  SCia< 

Y  ALS0*Y 

TP (ImCW.EG.OJ  GO  TO  5* 

GO  TO  68 
56  MPgsNR*'  +  1 

PPSOlMPfc,  1 )*X 1 
^PSY(^«'Ni,?)SYl 
I0O'IPL*0 
TFLAG35 
58  ISPi*I5 
ISaIS-1 

IMIS,FO,0)  ISsNS 
GO  TO  SO 

60  Nfi  Y*M«fvM 

RPSNliMGN.  1  )•*? 

BP 8N(NBM,2)«Y? 

IOOURLaO 
IPLAG«5 
68  ISsISPI 
1 3  P I  ■  I R  ♦  1 

IPUS.FQ.NS)  I5Pl*1 

K  0 1 P  a  1 
GO  TO  50 

70  IF(r«RPL3Ni,F0,O>  GO  TO  74 

CALL  CHECK2(X,Y,RP3Ir<‘,«i«Pl8N,NHPTSN,ICHK2) 

IP(lCHt<2,EG,n  GO  TO  120 

74  OPPLSNbNRPLSN+1 
GO  TO  76 

75  CALL  CHF.CK1  IX,Y,PPSM,MHN,ICHK) 
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IFnCHK.F-O.l)  GO  TO  90 
76  NRNa\B*+l 

BPSN(NRN, 1 )sX 
«PS'»(MRNf  2)SY 
X  1  sX 
V 1  a  Y 
IOOUBLal 
XAtSOsY 
YALSOsY 

IF  ( IFLAG,Fo,6)  TFLAGal 
IF(IFLAG,HJ,1  )  G H  TO  35 
IBPiaKLfl 

iFfKL.f Q.NR)  IPPJ=1 

TFLAGa/J 

GO  TO  IS 

40  IF  (FjRF',fc«,n)  GO  TO  1?0 

IF  (  ■4 W *m # i_ f  ,2)  *RPiS#«.vRPLSr*-l 
IF  (NRM.l F ,?)  GO  TO  1?0 
IF  (NEXT ,F 0,0)  GO  TO  100 
NRPTSNff'RPtSMsNRN 

CALL  TH  A^SP  ( WP8I0»  kWPLSN,  SP31  , NS , RP8N , SP3 , 3 ,GAMV , GAML  ) 

WRl'ISO 

x 1 sXSLO 

YtaYSFP 

I F  L  A  G  ■  6 

I RP 1  a  OF  XT 

lOOOMLel 

XALS0SX8F P 

YALSOaYSF P 

I F  ( ^WPl 3N , fcO , S  >  GO  TO  120 
GO  TO  XS 
100  x 1 s  S  P  S ( 1  ,  1  ) 

YlsSF’Sfl  ,2) 

CALL  CTRCLf  (kPSGJip.l,  xl  *  Y1  ,  I.YCR) 

IFiriCP.F  0,0)  GO  T i •  1  10 
X1SRP3  1  ) 

YlsRP3N(f  Rfv,?) 

CALL  CLOSE  (  SP  S»NS,X1  ,Y1,HM 
Xl=RP3*  (1,1 ) 

YlsRRSf  (1 ,2) 

CALL  Cl  OSF.  (  SP3,  nSf  x  1  *  Y 1  » *  )  ) 

Oil  103  kO I Rs  1 ,2 
K  a  K  0 

IF aCTsf-1 )**(K0IR+1 ) 

on  102  1=1, ns 

Lstf  +  IF  A  L I *( 1-1) 

IFfL.GT.fS)  Lai 
IF  fL.LT.I  )  L*NS 

I  F  ( L  , F.  0 , 1  t  Afi[;#KOIR,F.U,1  )  KC2*I 
IF  f  l  , F  0  ,  f'S ,  AX'f>, KO  I  P  ,  F  0, ?)  KsoS*l-1 
MsMRli+  I 
00  101  1=1 ,2 
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101  HP8N(*,J)*3PS(L#J) 

IF(L,EO,M)  GO  TO  1020 

102  CONTINUE 

1  020  00  1025  KKs 1 #  NS 

i  1 aSPS ( KK 1 1  ) 

V 1  *8PS ( KK ,  ? ) 

CALL  CIWCLf (RPSN,M,X1 , Y1 ,  INCR) 

If (INCR.FQ.l )  GO  TO  10S 
1025  CONTINUE 
GO  TO  10a 
10*  CONTINUE 

vHITF (6,200) 

200  FORMAT  (1U1  , 'SOMETHING  ^HfiNG  3HADOm  SUBROUT  I NE  *  ) 

STOP 

10U  NRP  TSN  ( IVRPL SN  )  sM 

CALL  I P  A^gP  ( RPS I  n  ,  nRPlSn  #  *■»  SP8I  #  MS#  RPSN  ,  SPS,  3 ,  G  A«V  ,  GAML  ) 
00  105  J=l,2 

105  RPS‘J(1  ,  J)*wPS(N«N,.n 
K*K'i 

DO  107  1*1, NS 
L*k-IEACT*(I-1  ) 

IFd-.KT.MS)  L*1 
If  a. I  T.n  L*NS 

If (L.EO. 1 .and, If  ACT ,E 0.-1  )  k*2-I 
IF(L,f>V  S.  AND.  If  ACT.f'-l,  1  )  KaN$*  I«1 
n*I  +  1 

DfJ  10h  .1*1,? 

106  PPSN(ui.n«8PS(L,J) 

IF(L,f'\M)  GO  TU  10H 

107  COOT  I NUF 

10«  m*^+1 

DO  109  J* 1 , 2 
109  *PS-J(m,  jjartPSO  ,J) 
v)9PLSO*:VWMLSN  +  1 
GO  TO  115 
j  j  o  ^snRm 

11  5  jPTsr  f  S'  )=M 

CAI.L  1  w  •.  SI  f"MM  '  ,  •  rtFl.3..,  ,SPSI,  S,«PSn,SPS,  l.GA^V.CA.vL) 

12  0  i*  •  •  ■  r  i  ■  ► 

n^i  •  t' i_  *■> 

Iff  >  ►  I  S  ,  I  ,  M  hi-  I  ■  IV 
•  -  J  M><  I.  p  1  ,  •  "  P|  S 
‘  K  S  0  I  S  ( I.  ) 

-r  o  1  s  (  I  )  =  - 
>.|1  MO  l  i  1  ,  -h 
nyswf.ST  U  ,1,1) 

YUSUPS T  f I , T , 2) 

ZusrfUST  ( L , I , 5 ) 

CALL  P|  '•  fl(i',SC,o,SI.,.S  ,  so,  U  )  SC  ,  *  U  ,  YU  ,  Za,  YX  ,  y  Y  ,  7Z) 

^psiil,  i  ,n=Y < 

RHSI  f L , 1 , 2 ) s V  Y 
Mo  RPSI(L,I,  J)*/Z 


WE  TURN  62 

E  on 

Cop 7  available  to  DDC  doO^f 
oermlt  fully  Icniblo  xcprodhtcMoQ 
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SUBROUTINE  INTERFCFPTS,0IRC08,NFENCE,NFPT8,Xl,Yt,Zl»REFLTN,ieLEAR) 

DIMENSION  APT (2? 5) , 8L (3) 

DIMENSION  EFTS  fS>  30737 ,  DT*COSf5,3) ,NPPT8t  55# 

1  Xl(3),XB<S>,REFLTNfS),COEPB<4),FPS<30,2) 

REAL»8  COEFBFC  5,4) 

COMMON /FENCE/ COE FBF 
DATA  TOL/1 ,0P*04/ 

REAL*#  0IRC08,X1,Y1,Z1 ,REFLTN,COEFB» SL,XB, XI,  X,Y,Z, DISC, DABS, DARSIN 

xrmixi 

XI (2)*Y1 
xnsjizt 
D0~s  m,r 

XB(I)*0,0 

SL<I)*0.0 

oo  r  jii.i 

XB(I)*XB(I)*DIRC03(I,J)*XI(J} 

5  8L(I)*8L(I)tDlRC08CI,J)*REFLTN(J) 

DO  50  FTH  ,NF£NCE 
DO  10  J«l,4 

10  C0EF8CJ)*C0EFSFCK,J) 

CALL  PL»NPTfCOEFBCl),COEFB(2),COEFB(3),COEF8(4),8LCl),8L(2), 

1  8L(3),DI8C,XB(l),XB(2),XB(3),X,Y,Z) 

1F(OABS(OI8C),LT,TOL)  GO  TO  50 

C»U  80tN0Tt8L-CnYSLt«)*St(3),XB(l),XB(2),XB(3),X,Y,Z,INT) 

IFUNT.EO.O)  GO  TO  50 

IP(INT,EQ,2)  GO  TO  55 

DI8T»0 , 0 

DO  12  1*1,3 

12  DlST*D13nC0EFB(I)*C0EPB(l) 

DI3T«SGRT  C0I8T) 

DO  14  1*1,3 

14  C0EFB(I)*C0EFB(I)/DI8T 

IF (DAB8CC0IFB (1 ) ) ,GT , TOL ,OR,OA03(COEFB (2) ) »GT , TOL)  GO  TO  16 
GV*i, 570796 
Gl«0,0 
GO  TO  !B 

16  GV*DAR$IN(COEFB (3) *0.9999) 

Gl*DARSIn(C0EFB(2) /COS (GV) *0,9999) 

IB  SNGV*3IN(GV) 

CSGVaCOS(GV) 

8NGL*SINCGD 

C8GL*C0S(GL) 

AFTCl,l)*3NGV 

AFT ( 1 ,2)*0,0 

AFT ( 1 , 3)*CSGV 

AFT (2 i 1 )*C8GV*SNGL 

AFT (2,2)*CSGL 

AFT(2,3)*"8NGV*8NQL 

NF*NFPT8(K) 

00  90  IF«1,NF 
DO  90  J*l,2 
FP 8(IF,J)*0,0 
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no  a  o  *  h  =  t  #  5 

an  FPSnF,J)»F**S{IF,J)*AF  T(J*KK)*FPT5(Kf  IF,KK) 
jMPaAFT(l.n*X*AFTCt»S)*Z 
YtP*AfT(?*l)*XtAET(?,<?)*YfAFT(2*3)*Z 
CALL  CI»CLK(FP9,wF,XlP#ViP,  I*mC«) 
if  ( l)  tn  Tn  S5 
SO  COoTlMt'E 
(?ETij««y 
ss  ICLfcAHel 
PE  TUP* 

F  .m  n 
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SUBROUTINE  QU ADPT ( A , B , C , D , E , F , G , H , CK , Cl , Si , SK , 8N , X i , Y 1 , Z 1 , X 1 1 , Y 1 1 , 

1  ZI1#XI2,YI2,ZI2#INTC) 

IMPLICIT  REALMS  (A-H,0-Z) 

3QRT(XX)«0SQRT(XX) 

A8S(XX)«0ABS(XXl 
DATA  TOL/1.0E-08/ 

INTCpO 

IF(ABS(SU,LT,ABS(SN).AND,ABS(SM),LE,ABS(SNn  GO  TO  30 
IF(ABS(SL),LT,ABS(9M) ,AN0,ABS(SN),LT.A8S(SM))  GO  TO  BO 
3A*A+B*(3M/gl)*#2+C*(SN/sL)**2+D*BM/Sl+F*SM*SN/SL**2+E*SN/SL 
SB«2,*B*Vt*S^/Sl-2#*8*Xl* (SM/SL)**2T2,*C*Zl*SN/SL»2,*C*Xl* 

1  (8N/8l)*A24D*Yl«0*Xl*8M/8L*E*Zl-t*Xl*8N/SLtF*Zl*3M/SL*F*Xl*SM* 

2  SN/SL«*2*F*Y1*SN/SL*F*X1*SM*SN/SL**2*G«>H*3M/SL*CK*8N/Sl 
SC«B*Yl*Yl*B*Xl*Xl*CSM/8L)**2TC*Zi*Zl-2,*B*Yl *X1*SM/SL+C*X1*X1* 

1  CSN/8l)**2-2,*C*X1*Z1*8N/9LyF*Y1*Z1-F*X1*YUSN/SL»F*X1*Z1*SM/SL4 

2  F*Xl*Xl*Sw*SN/Sl**2*H*Yl»H*Xi*3M/SL+CK*Zl-CK*Xl*3N/SL*CL 
ICONTRil 

§  DlSC»SR*SH*-«, *SA*SC 

C  WRITE (Nf 200)  ICONTR,  X1,Y1,Z1*SA, 

C  1  38, SC, DISC 

C200  FO»MATnM0,I5,5X,6Et5,«/lHo,10X,«EI5ta) 

IF(0I8C.lTt0AB8(SA)/10,)  RETURN 
IFfOISC.LT.O.O)  RETURN 
IF(AHS(3A).GE.T0l)  GO  TO  10 
1F(AB8(SB),LT.T0L)  RETURN 
X 1 1 *»SC/SB 
I N  T  C  » I 

IF C  ICONTR, FO, ?.)  GO  TO  *0 
IF(ICONTW,EQ,3)  r,0  TO  «0 
GO  TO  ?0 

1 0  Xll*-8B/8A/2,*S(JRTCDI8C)/SA/2, 

X 1 2»"8B/SA/2,»8QRT  (OJsO/SA/2, 

INTC«2 

IF ( ICONTR , fO , 2 )  GO  TO  bO 
!F(IC0*TR,FQ,3)  GO  TO  RO 
20  Yll*YlTSM*(Xl1»X| )/8L 
ZIl«ZU8N*(XIl.Xn/SL 
IF(tNTC.FD.l)  RETURN 
YI2«Y1+SM*(XI2*X1 )/8l 
ZI?«ZI+8N*(XI2.X1)/SL 
RETURN 

30  S*»C +B*  (SM/SN)**2TA*(St/SN)  **2  +  F *SM/SN+f  *3L/sN  +  l)*SM*8L/SN**2 
SB«2,*P*Y1*SM/SN-2,*B*Z1 * (SM/8N)**2*2,*A*Xi*SL/SN«2,*A*Zl* 

1  (Sl/SN)**2  +  F#Yl-F*Zl*8F'/SN+E*Xl-E*Zl*Sl/8NT0*Xli»SM/SN-0i»Zl*SM* 

2  8L/SN**2  +  0*Y1  *8l/8N»0*Zl *SM*SL/8N* *2  +  Ck +H*8M/SN+G*Sl/SN 
SC3«*Y1 *Yl tH*Zl *Z1 * (S^/Sn) **2-2.*B*Yl *Z1  *SM/8N  +  A*xl *X1 +A*71 *Z1  * 

1  (Sl./SN)**2-2,*A*Xl*Zl*Sl/SNTD*Y1*X!*n*Yl*Zl*SL/SN-D*Xl*Z1*SM/SN* 

?  0*Z1 *Z1 *Sw*SL/8Na*2*h*Y1 «h*Z1 aBM/SNTG*X1»G*ZI  *SL/SN*Cl 
ICONTR* 3 
GO  TO  5 
ao  Z 1 1 » X  T I 

YI13Y1+S^(ZI1-Z1  )/SN 
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Xll  ■X1+3L*(ZI  1-ZD/3N 
IFdMTC.EQ.l)  RETURN 
ZI2»XI2 

YX2«Y1*SW*(Z!2-Zn/SN 

XI2«X1*SL*(ZI2-Z1)/SN 

RETURN 

50  SAaB+A*(3L/SM)»*24C*(SN/8M)**2f0#3L/3MfE*SL*3N/3M**24F*8N/8M 

S8B2t*A*Xl*SL/SM*2,*A*Yl*(SL/SM)**242(*C*Zi*SN/S^o2.*C*Yl* 

1  (8N/3M)**2«D*Xl«D*Yl*3l/3M*E*Xi*3N/3M*E*Zt*3L/3M-2,*E*Yl*3L*3N/ 

2  SM**2+F*Zl«F*Yl*3N/8M*6*SL/3MfH+CK*8N/SM 
SC»A*Xl*Xl+A*Yl*Yl*(SL/SM)**2",2i*A*Yl*Xl*SL/8M+C*Zl*Zl+C*Yl*Yl* 

1  (3N/3M) **2*2, *C*Y| *Zl *3N/3^YE*X| *Zi»E*X| *Y| *8N/3M»E*Y1 *Zl *3L/3M» 

2  E*Y1*Y1*SL«3N/3M**2«G*X1«G*Y1*3L/SM«CK*Z1«CK*Y!*SN/SM*CL 
IC0NTR«2 

GO  TO  5 
60  Y 1 1 *X 1 1 

XI1«X1+SL*(YH-Y1)/S* 

Zl 1 »Z1 +3N*  t Yl t  »Y1 )/SM 
IF(INTC,EO,l)  RETURN 
YI2«XI2 

XI2»X1*3L*(YI2-Y1)/3M 

Zl2«Zl*3N*(yI2-Yl)/3M 

RETURN 

ENO 
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SURHOMTP'F  SnRNnT(SL#SM,SN  .Xt  ,  Yl#Zl  »x*y* 
IMPLICIT  W  fe  A  t..  *  R  (A«H,D"Z) 
ARS(KXlanARS(YX) 

OATA  TOL/1  .Of’  -0A/ 


IShAHs? 

SU*X-Xl 

3Miav«vl 


IF(UBS(SLnTARS(Sm)tAHS(SNm.LT.T()U 


I F  (  ABS  ( SL )  ,GT  ,TUU  DIR  f l_l' -L  /SM 
TF(AHS(S^I.GT.TnLl  0 1  MSl/ T  R  +  ^ 
IF(A«S<3M,GT,T0l)  0IR-0IH+3N1/3N 


I Sh AOs  1 

IF (DIR,  IT. 0.0)  ISHAOsO 


RETURN 

END 


2, IShAD) 


RETURN 


St'RWOUT I  ME  PR0HL(ALPHA2,RFT A2,GAMMA,  NPSI ,NDISP, NRPTS, 

1  Y9CMIN,  YSCMAX,8UMk,  IGAM,NPANEL) 

OImF.NSION  ALPHA?  n  0001,  HfcTA2(  1000)  , NRPTS (7, U,  t 8 ) , JmNMX ( 1 0 , 2 ) , 

1  HPTSU,?),8PSC50,?),PTS(30,2) 

COMMQN/MLPH/vGAMS(7,a) 

AHEAfXl , Y1 ,X?,Y2,X3, VS) =ABS ( ( X J * Y2+X2* Y3*X 3* Y 1 «Xi* Y2» 
l  Xl*Y3«X2*Yl)/2,) 

TF (NPSI ,FQ, 1 )  RETURN 

YRTTE (8,200)  GAMMA , YSCMJu, YSCMAX 

200  FORMAT ( l hi  ,  o  9  x , 'PROBABILITY  EUnCTIUNS'/1H0,«3X, 'RATIOS  OF  SWEPT  AR 
1 E  A  TO  TOTAL  ARF A '/ 1  HO , bX , 'SUN  ELE  V  A  T  I  ON  '  ,  F20 , 2/ 1  H  ,5X  ,  'MINIMUM  OSS 
2ERVER  ELEVATION', F7.2/1H  ,bX, 'MAXIMUM  OBSERVER  ELEVATION', F7.2// 

3  |HO,iOX, •SURFACE  NO, ' , 5X , 'SWEPT  AREA', SX, 'PROBABILITY  FUNCTION') 
00  2000  KL«1»NPANFL 
TOTAaO.O 

M*MGAVS(KL, IGAN) 
no  20  L *  1 ,NPSI 
NPT3«NRPTS(KI.,IGAM,L) 

TMNMxfL , 2 ) *0 
JJsO 

00  20  k*1,N0ISP 
00  20  via  1  , nPTS 
m*m*1 

IF (ALPha?(M) ,GT,90, )  GU  TO  20 
JJ«JJ+1 

IEUJ.GT, !  )  GO  TO  10 
IFCL.GT, 1  )  GO  TO  b 
RminsHE  TA2(m) 

HMAX*HMIN 
5  IMNMX(L,1)»M 

IMNMXIL  #2)«M 

YMAXiAI.PHA2(M) 

YMINsYMAX 
GU  TO  20 

10  IF(ALPHA2(M),LT,YMIi\i)  IMNmX(L, 1)»m 
IF(ALPhA2(m),c.T,YMaX)  IMmmx(L,2)«m 
IF(MtEO,IMNMX(L,l ) )  YMlNaALPHA2(M) 

IF ( IMnmX (L»2) ,EG,M)  YMAXaALPHA2(M) 

IF(8F.TA2(M)  .LT.BmIN)  RMINbhetA2(m) 

1F(0FTA2(M) .GT.BMAX)  RMAX«BETA2(M) 

20  CONTINUE 

sPsn,n»BMiN*5, 

SPSd  ,2)«YSCMIN 
SPS(2, 1 )«BMAX+5, 

SPS(2,2)»YSCMIN 
SPS(3, 1 )«flMAX*5. 

SPS (3,2 ) a YSCMAX 
SPS(«,1)«HMIN«5, 

3PS(«,2)«YSCMAX 
NMlaNP8I*l 
no  looo  l*i,nmi 

IF(IMNMX(L,2)tEQ,0)  GO  TO  1000 
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■"  >V’**rt'  t^vn**  .►.  v-r!*-  *  r*t*r».  -*w\fv  'i+c 


LPl*LM 

IF (IMNMX  CLPt »?) ,EQ,0)  GO  TO  1000 

no  30  1*1,? 

MaTMNMxCLf I) 

BPTSd  ,  1)«PKTA?(M) 

30  BPTS(I, 2)«ALPHA?fM) 

no  ao  i*i,? 

N  *  3  •  I 

^■TMNMX  ((.PI  ,  I  ) 

BPTSI-M  +  ?,  1  )aPFTA?(M) 

UO  BPTS(N  +  ?,?)*A|.PHA?(M) 

IDOUPL*0 
ICHK*o 
LH'jtso 
m  P  T  9  *  0 

I«*t 

50  XCsOPT  S ( I H,  1 ) 

VC*BPTS( IW,2) 

CALL  CIWCLf (9PS,a,XC, YC, IMCM) 

TFCINCP.PE,?)  C,D  TO  HO 
IF(I»tEO,«)  GO  TP  60 
IPsTP+1 
GO  TO  50 

h o  no  70  l*i, a 
no  70  jsl ,2 
70  PT9(I.J)*«PT5(I#J) 

OPT  Sea 
GO  TP  900 
MO  I FL  &Gs i 

IF(IM>',feU,0)  IFL  AG*2 
IF(IFLAG,t«e2)  GO  TO  MS 
MPT8*1 

PT9(WPT8,  n*BPT8dH,  1  1 
PTS(^PTS,2)sapTSf IP,?) 

*5  IWPlsIw-M 

IF(Ift,k«,a)  IWPlsl 
90  XlsBPTSCI»,l) 

YlsMPTSUW,?) 

X?eHPT  8 ( 1  BP  1 ,  1  ) 

Y2sBPT9( I»P| ,2) 

CALL  I  I'TF  wC  (  SPS,  a  ,  X  1  ,  V  1  ,  X?  ,  Y?  ,  X  ,  Y  ,  -JL  ,  IF  AULT  ,  muUHl. , 
1  XALSO,  YALSCl) 

IF  (  U  A  *  1 1.  T  ,feQ,0)  GO  TO  110 

IF(IFLAG,FtJ,2,OH,lFLAG,fcri,3)  GO  TO  1 00 

IF  f  •'!  P  T  ? , UT , 0 )  CALL  ChFCKl  (X2,Y?,PTS,NPTS,  ICHK) 

IFflCHK  ,Ffj,l )  GO  TP  900 

•MPTSsoPl  9+1 

PTS(NPTS, 1 )*X2 

PT9£^PTS,?)*Y? 

IDOilMLrO 
t Ft  AGs  1 

100  IF  dFLAG.Fo,  31  IFLAGs? 
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JS*l*‘**PW*<Tr  '  "  -  wwtravw 


IFCIRPl.EQ.i)  CO  TO  900 

105  I R» I R ♦ 1 

IRP1«IR*1 

IF(IR,EG.«)  IRPl«l 
GO  TO  90 

110  IFHPTS.GT.O)  CALL  CHECK  1  ( X ,  Y ,  PT  3,  NPTS#  ICHK ) 
IF ( ICHK ,EQ, 1 )  GO  TO  900 
NPTS«NPT8M 
PTS(NPTS,1)«X 
IF(IFLAG,EQ,«)  IFLAG*2 
IF(IFLAG,EU,2)  GO  TO  105 
IOOUBL*1 
XALSOaX 
YAl30«Y 

IFCIFLAG.EG.I)  IFL AG«3 
IF(IFLAG,E^,2)  IFLAGaA 
GO  TO  1  OS 

900  IF(NPTS,LE,2)  GO  TO  1000 
N*MPT8«2 
X 1  aPTS ( 1 » 1 ) 

Y  1  aPTS ( 1  #  2 ) 

00  950  Isl,M 
X2*PTS(  I  +  l # 1 ) 

Y2»PT8(l4t,2) 

X3*PTS  ( I +2i 1 ) 

Y3*PTS  T 1 *2»  2) 

950  TOTA«TnTA*ARfcA(Xl ,Y1  ,X2,Y2,X3,Ti) 

1000  COMTI'MUt 

PRTOTaTOTA/SUMfc 
IF (PWTOT.GT, 1 ,0J  PRTUT«1,0 
WRITEC6.201)  Kt.  »  TOTA ,  PRTOT 

201  FO«MAT(lH0,3«X.'(',Il,')'.9XrF7,2,16X#F5,3) 

2000  COMTINIJE 
RETURN 
E  NO 
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SUBROUTINE  PLOT (XPTS,YPTS, GAMMA, IPANEL»NPSI,NDI3P,NRPTS, 

1  X3CMIN,X3CMAX#Y3CMIN#V3CMAX,PHI,THETA#1GAM,NPANEI#IANGI) 

DIMENSION  XPTS(  1  000),VPTS(1000)#NRPT3C7,4»,18),3YMBOLU8) 

DIMENSION  XPH(10),OUT(91) 

DATA  SYMBOL/ 1  HI ,1H2,1H3, 1H4,JH5» i H6, 1 H7 , 1H8 , 1H9, 1 hA,  1 HB, 1 HC , t HD, 
i  IHE,1HF,1HG,1HH,1HK/ 

DATA  BLANK, PERI0D#FI/1H  ,1H,,|HI/ 

COMmON/PLPH/MGAM8(7#«) 

IFdPANEL.EQ.NPANELl  WRITF(6,?08)  IPANEL 
IFdPANEL.NE.NPANEL)  J TF (6 , 20* ) 

WRITE(8#200)  GAMMA, THETA#PHI 

DELX«(X3CMAX*XSCMIN)/R0, 

0ELY«(YSCMAX-Y3CMlM/«0, 

KPRINT«S 

TY2aYSCN4X+DFLY/?, 

DO  1«0  1*1, 41 
TY1*TY2»0ELY 
DO  SO  J«1 ,91 
50  OUTCJ)*BLANK 

IF(KPRINT,F.0,5)  nUT(l)*FI 
IF (KPRINT.NF.5)  0UT(1)*PF«I0D 
IF(*PRInT,EO,5)  OUT (91 )*F I 
IF(KPHInT.ne,5J  0UT(9t)*PFRI0D 
IFd.NEtl.AND,I.NE,«n  GO  TU  50 
L*9 

DO  UO  J* 1 ,91 
L*L  + 1 

OUT (J)sPF  HIOD 
IF(L,FN,lO)  OUT ( J ) c F I 
40  IF(L,EO,10)  l,*0 

50  00  90  KL*IPANEL,NPANFL  ‘ 

M*MCflMS(KL, IGAM) 

DO  90  L*1  ,NPvSI 

nPTS*NRPT3(KL,IGAM,L) 

DO  90  K*1 ,NDISP 
DO  90  J*1,NPTS 

M*M*  1 

IF (VPTSCM) ,GT,90. )  GO  TU  90 

Y*YPT3(M) 

IFLAGsl 

IF(J,E0,1)  ^Ml*MfNpT3*l 

MPl *M+t 

IF(JtEQ,NHTS)  mpj*m.npTS+1 

TF(I,E«.1  .AND,v,r,Ft(V8CMAX+f)ELY/2,))  GO  TO  52 
IF(I,EO,ai  ,Ai'0,Y,LTt(YSCMlN.DELY/2,))  GO  TO  56 
IF(Y,GE , TY2.0P, Y.LT.TYl }  GO  TO  90 
xaxPTS(M) 

GO  TO  BO 

52  IF(YPTS(MMi),GF_,(YSCMAX*DfcLY/2.))  GO  TO  54 
IFlAG*? 

X*XPTS(*V1  )  +  (ySCMAX-YPT3(NMi )  )*(XPT8(M)*XPTS(MM1  )  )/(YPT8(M>- 
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1  YPTS(  v'Ml  )  ) 
go  ro  *o 

Stt  IP  ( VP  T  S ( VP 1  )  ,  GP  ,  (  YSC  *'fl X ♦('  fe L  Y  /  a ,  )  )  r,n  TO  RO 
IP  l .  A G s  S 

X  s  X  P  ?  S  f  v<  P  j  )  ♦  (  YSC*  ax-yPTS(NP1  )  )  *(  xPTS(r*)«XPT8(MPi  )  )/(  ypTS(P)- 
i  vprs(“pn) 

GO  TO  NO 

SN  TF(YPT3('«nn.LF#(Y8C'‘'TM-nFLY/2,))  GO  TO  SB 
IFLAGa/4 

X3  XPT  S  ( ”  )  M  vSCf'  I N- YPTS  ( * )  )  *  (*P  TS  (  “*  '  t)  -/P  T  S  M) )  /  (  YPT  S  (MV  1)  * 

1  YPTST  )) 

GO  TO  M) 

SB  IP  (  VPTS(f'Pn  ,UP  ,  (YSC'l^-OM.  V/2.J)  GO  TO  VO 
l  fi  agss 

XaXPTSf *)♦ (YSCVTN..YPTS(N) ) *(XPTS(PP1 )-XPTS(M))/(YPTS(MPl)- 
l  vp  IS  (*  ) ) 

Ml  I  X?svsr  MlMfDH  */? >, 

00  70  k  K  s  t  ,  Q  1 
T  X  J  aT  X?»DF|.X 

IPfX.GP tTXl ,AMn#x,LT,TK2)  GO  TO  B 0 
70  TXPsTxPfOflX 
BO  LI  =**  +  ( 1 -I AUGl ) *  f ) 

OuT(L.UsSV-.ROL(L) 

GO  Tn  (90, Su, 90, SB, 90),  I P  L  A  U 
90  CONTINUE 

VPLOTsTY?-r'El  Y/2, 

rP  <KPPTMtFQ,5)  ■vRTTfe(*#?01  )  YPLOT  ,  OUT  ,  YPLr.T 
TP  (*PRTf  T.'  fc  ,S1  -«lTf(6#?i)S)  OUT 
IF  (*Pr"T>«TtlrGtS)  kFHMsH 
K  P  R I f  •  T  e  K  P  W  I N  T  ♦  1 
lao  rY?sTV?-nELY 
XPP(  1  JaXSC1'!  M 

IF  (I  AMGL.FG.O)  XPH  ( 1  )sX8CMAX 
1)0  1  SO  T  r2  ,10 

ISO  XP9(  I  )SXPR(  F-1  )  +  n*T  A'^Gl  )  *1  0,  *  OF.  LX 

^  w I T  p (b'jau) 

VWITF(F,?0?I  XPR 

a«ITE  (N.P07) 

*06  FOPPATClH|,iinX,'SUP  GLINT  SIUPATimF  FOR  SURFACE  NO,*, IS) 

?00  FQ^iTflPfl,  *Sl'N  f  l  P  vat  ion*  , 

1  £12,2/1*  ,  'A/C  PI  TCH  ATI  I TUOF *,F7,2/IH  , 

?  *A/C  POLL  ATTITUDE'.  FM.2//1H  .'OBSERVER  ANGLE '  x  93X , 
h  *ORS£PVFP  A  * GLF  *  / 1  *  , 3X, 'DEGREES*, 100X,* DEG WEES*/) 

?07  FORMAT  ( / »  SR  X ,  *GLINT  AZIMUTH  *IT0  RESPECT  TO  SU.v  -  DtGRf.ES*) 

P01  F 0»m A T ( 1  m  ,F10#5,5x,91A1,1 x,F10,S) 

aos  P  ORMAT  ( t  M  ,lSX,9jAn 

?0?  FORMATHH  ,9X,tOF1ft,S) 

aOU  FORMAT (/) 

209  rorm  A  T  ( J  Hj ,  97  X ,  *TOTAL  SO  N  GLINT  SJGNATUKFM 
RETURN 
Eno 
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SUHwniiTP'E  Tt»AMsP(»Psl#IkPL,NH,8F8l#N8#»P8#SPSf  IC#GA«V,OAML) 
ni^f-fiSTOv  WPS  I  (S.  50,51  ,$PSI (30* 31 ,RPS( 50  * ?) * 3PS( 30 ,2) 
COwMf'M/  TWA^'/AT  f  2 »  5 1  #ATff'V(5,2) 

r,n  r  n  no,  tin,  tin),  ir. 

1ft  $;-irivsSTf'(t;rtMv/S7,3) 

f.  sc;  \/  =  r  n s  ( r,  n  m  v /«s 7 . 5 ) 

Si'iRL  =  SIiMCr,Ar..i  /S7.3) 

CSRL»CnS(r,AFiL/57t  5) 

at  ( i ,  i  )sS  vh;v 

ATM  ,?)=0,0 

AIM.  51SC3GV 

AT  (?,  1  )aCSGV*3i'<Rl 

AT(2,?)sCSG|. 

AT  (2, 3)«*SA!fiV*SNGL 
•  >n  is  T»1 » 3 
t»0  IS  JM  ,? 

IS  A  T  !  M  V  (  T  ,  J  )  s  A  T  (  T »  T  ) 

•  OU  tOO  I r  1  ,  >*>i S 
00  1O0  Jsl ,2 
SPSf I, J)*0,0 
00  too  K*1#3 

100  SPS(  I  ,.f)sSPS(  T,  JHAT(J,K)*SPST  (I,K) 

110  00  120  Ist,N.-« 

no  120  J*1 .? 

WPS  (I  , J)«o,n 

00  120  KS1 .3 

120  PPS(J».1)>RPS(I»J1«AT(J*K)*PPSI  (  I W  P  l  ,  I  *  K  ) 

RE  TilWW 

1  JO  00  100  I  a  1,  IMP 
00  1  an  .1  ■  l  ,  3 
WPS  T ( I  PPL . I , J ) a o , 0 
00  tan  K3t,2 

iao  WPS  I  ( rwPl. ,  T ,  J)«RPST(TRPl,I,  J)*ATlMV(J,K)*PPS(Ii  K) 

RET.iWM 

Mil 
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SUBROUTINE  CH£CK?m  #  V 1 # ARRAY # NPL8# NPTAR# ICOMP) 
DIMENSION  ARRAY(5»30#3)  #NPTAR(5) 

COMMON/ TRAm/ AT  C2#3)#ATINV(3#2) 

OATA  TOL/1 ,OE-Oa/ 

ICOMPao 

X2»ATlNV(l#D*)fl+ATINV(l,2)*Yl 
Y2*A  T I M/  f  2#  l)*Xt  4ATINV <2# 2)* Y! 
Z2«ATINV(3,nAXl+ATINV:J,2)*Yl 
00  10  LSI ,NPl s 
mPTS«NPTAR(L) 

no  10  Ial^'PTS 
XaARRAY (L# I  #  1  ) 

YaARRAY (L  # I #  2) 

23ARRAY (L  # I #  3) 

IF (AflSfX2«X),GE,T0L)  GO  TO  10 
IF(AH8(Y2-V)#GE,TOl.)  GO  JO  10 
1FCABSCZ2-Z).GE,T0L)  GO  TO  10 
I COSP« 1 
RETURN 
10  CONTINUE 
RF  T'lWN 
F  NO 
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SUBROUTINE  CLOSE C ARRAY, NPTS, XI ,Y 1,10 
DIMENSION  ARRAY ( SO , 2) 

K*1 

DIST*( ARRAY (1 , 1)-X1)*»2*( ARRAY (1,2) -Y1  )**2 
IFCNPTS.FQ.D  RETURN 
DO  10  I«2,NPTS 

oist2*(array(i,  n«xn**2f(  array  <i,2)  »yi  )**? 

IF (0IST2.LT .DIST)  K*I 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  INTE*C(P0INTS,NPTS,X1,Y1 #X2,Y2,X,Y,KL#IFAULT, 

1  I00UHL,XALS0,YAL30) 

DATA  TOt/1 ,0E-04/ 

DIMENSION  POINTS (30 *2) *  I PTS ( 30 ) # XJnT (30 ) , Y I  NT  130) 

NlNtC*0 
IF  AULT* 1 
00  7ti  1*1  *  NPTS 
IPt*I*l 

IF(I,EO.NPTS)  IP1 ■ I 
X3*P0INTS(I,1 ) 

Y3*P0InTS(I»2) 

XAbPOInTSC IP1 i 1 ) 

YAsPOINTSCIPl ,?.) 

XM  J  N*  X  3 
XMAXBX3 
V  M  [  'v  8  ¥  J 
YMAX*Y3 

TF(X«,LT,X5)  XMIN*X4 
IF ( XA  t GT, X3 )  XMAX*X« 

IFCYU.LT.YJ)  YMINbYU 
IF(Y<i,GT,Y3)  YmaxbYU 
0EL1*X2-X1 
OEL2*X«-X3 

IF(A0S(DELi),LT,TOU  GO  TO  10 
S A 1  * ( Y2*Y 1 ) /OEL 1 
SP1*Y1-SA1*X1 

10  IF(A0S(OFL2).LT,TOL)  GO  TO  20 
8A2S(Y<UY3)/U£12 
S02*¥  S-SA?*X3 

20  IF(APSrOELn.GE.TOL)  GO  TO  AO 
IF(ABS(OEL2)§Ge,TI)L)  GO  TO  30 
GO  TO  7 A 

30  NlNTCoNINTCtt 
XIMT(NINTC)*X1 
YInT(NInTC)«8A2*X1*8»2 
IPTS(NTNTC)»I 
GO  TO  7n 

ii 0  IF(ARS(l)EL2),GE,TOL)  GO  TO  SO 
N I  M  T  C  *  N I  N  T  C  ♦  1 
X  I  M  T  (  N I N  T  C  )  *  X  3 
Y I N  T (NI VTC )*SA 1 *x3  +  SP  1 
IPTS(NINTC)*! 

GO  TO  70 

50  iJIFFA*8A1-SA2 

!F(AHS(OTFFA),LTtTOl.)  GO  TO  7<J 
NlNTC*OINTC*l 

XlwTfNlNTC)*(SP?-SH1 ) /O  IFF  A 
¥  I N  T ( N I  NTC ) *SA 1  *  X I  NT ( N I  NT C ) +SB 1 
IPTS(NIMTC)«I 
in  CONTINUE 

CAl  L  CMECK3{XMIN#X^AX#  VF"IN,YfAAX»XlNT(NINTC)  #  YINT(NINTC)»IPETW) 
IF(lHF.T«,EG,n  IPTS(N1NTC)*0 

76 


I 


fjt^VnQiu*. 


r 

k 


f 

j 

l 

i 

i 

{• 

l 
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74  criMT  I  fH.iE 
IF(WIMTC,GT,01  IF  AIJLTaO 
IF  (IF  AUI  T,F«,t)  PFTUPM 
X^IrjsXl 

x  N'  a  x  s  x  i 
Y^TNsYl 
YVAXaYJ 

IMXP.LT.X1)  XMj^ax? 
tFfX2.GT.Xn  X  M  A  X  *  X  2 
IF(Y?,LT,Yl)  Y  f"  I  f-.  =  Y  <? 

TFfYp.GT.Vl)  YOAXSY2 

no  50  T*1,A'IMTC 

if ( ir>m***t  ,fcQ,<n  r,o  to  75 

IFCAHSf  X  I m  T  ( !  5-XALSO)  ,IT,  TOL,AfMD,AHS(YlMT(I  )*YALSO)  ,LT. 
1  TOL.)  IPTSCTJ80 

75  IF(XT»»T(1),l.TtX*JN)  !PTS(I)aO 
IF  (XTMT  (I)  ,GT,X*AX)  IPTSdlsO 
IFlYlNtm.l  T.YMTU)  IPTS(I)  =  0 
IF(YTNTmtr,TtYMAX)  !PTS(I)bO 

fto  coMTiMoe. 

IFAULTal 

f>0  90  7  =  t  .NINTC 

IKU’TSm.FO.O)  GO  TO  90 

tF (IFAMLT.fcQ.O)  GO  TO  *5 

IFAOLTrO 

OISTs(Xp.T(T)-Xl  )***♦<  VIM  fI)-Yl)  **?. 

>*■1 

r,n  TO  9  0 

«5  0IST2afXlMT(I)-X1)**?+(VIM(i)-Tn**a 
IF(l)IST2,LT,0I8n  *«I 
IF(M,fcn,I)  0ISTSOI5T2 
90  C0MTINOF 

IFCIFAULT.fO.J)  RFTIIPM 
Kl»tPT9(M) 

X  a  x  I  w  T  ( M ) 

Y*V  JM  f M) 

RFTUPM 

Ffvl) 


77 


S'MriOiJTI.vf  OtTK-i(*MT!\.#XvAX,  fMlN#  Y*AX,X,  i,  I»ET*V) 

T  ><F  r- ■sr. 

IF  f  <,LT,yMtfO  IMF  T*  s  1 
IF  (X.riT.X  TflFT.»  =  t 

If  fY,!  T.YMM  T W f  T 5  1 
IF f Y.GT,  YVfiX)  MfcT/sl 
rfE  T'JWh 

H-y.i 

i 

i 


i 

i 

i 

i 
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SUBROUTINE  CIRClE(P0INT3,NPT3,XC# YC,INCR) 

OIMENSION  POINTS(30,2) 

OATA  TPL1/1,0F*0«/ 

NPC*0 

NNCaO 

00  10  la  1 # MPTS 
IP1 «!♦! 

iF(I.tn.NPTS)  !P1 ■! 

0I3Xl*Pt)INTSCI,l)-XC 
0ISY1sP0JNT3(I,21-YC 
OlSX?aPOTwTSC IP1,1)*XC 
DI9Y2*P0INTS(IP1,?)«YC 

IF(ABS(UlSXnfLT,T0L1  .AND, ABSCOISYIJ.LT, T0L1)  GO  TO  20 
TF(A«S(DlSX?),LTfT0Ll.ANn,ABS(DlSY2),LT,T0Ln  GO  TO  20 
ifc ABsroisxn.GE.Tnm  go  to  3 
IF(APS(OISX?) .GE.T0L1 )  GO  TO  1 
IF ((DISY1/OISY2).GT,0.0)  Go  TO  iO 
GO  TO  20 

1  IF(ABS(i)TSY2)  .GE.TOLl  )  GO  TO  2 
X*IHSX2 

GO  TO  p 

2  IP  (ABSfi)ISYl-OISY2),LT,TOLl)  GO  TO  JO 
GO  TO  5 

3  IF(AP8(UI8Yl),GE,TOLl)  GO  TO  32 

IF(ABSfOt8X2),LT,TnLl)  Go  TO  10 

IFfAP8(i)I3Y2)tGF1TOLl)  GO  TO  10 

IFC(f1ISXl/DISX2),GT,0,0)  GO  TO  10 
GO  TO  20 

32  IF (APS(ni8X2),GE.TOLl )  GO  TO  3« 

GO  TO  3h 

3 a  IF(AR8(OTSY2) .GE.TOLI)  GO  TO  36 
X  o  f )  I S  X  2 
GO  TO  « 

36  rF(ABS(OISY2-niSYntLT.TOLl)  GO  TO  10 

TP  (Ar38(i>TSX2*0lSXi  )  .GF.TOLl)  GO  TO  5 

IF((0ISV1/DI8Y?),GT,0#0)  GO  TO  10 
XaOlSX? 

GO  TP  8 

S  SAsroiSY?*Dl9Yl)/(OlSX2»f>l8Xl ) 

9Hat)lSYl«SA*OISXl 

Xa-SB/SA 

TP (X.LT.OISXi.AND.x.LT, DISX2)  GO  TO  10 
IFIX.GT.OISXi.ANP.X.GT.OISX?)  GO  TO  10 
0  IP  (ABSfXJ.lT.TOLl)  GO  TO  20 

IP  (X, IT, 0,0)  MNC.sMNC  +  1 
IF(X,GT,0,0)  mPCsnPC+1 
10  CONTINUE 

INCWsO 

IP  f 2* (NPf/2) ,E0,nPC )  RETURN 
lF(2*(MNC/2)iFQ,^NC)  RETURN 
IMCRsl 
RETURN 


20  INCH*?  79 

RETURN 
END 


SURtfrijTINE  CHECK!  (XI  ,  VI  ,  ARKAY,  ICOMP) 

DIMENSION  APPAY(S0,2) 

DATA  Tm/1,0^-0  U/ 

ICn-tp*r 
no  1ft  T  a  1  #  i'1 
*3 APPAY (  1,1) 

XSAPPAy (  T  ,?) 

IF(AHS(X.Xl),LT,TO|.f  A  NO , AHS ( Y. Y 1 ) , IT , TOl )  GO  TO  IS 

CO"'  r  timhe 
PF  ]  Ijk;-| 

Tfn-'Pst 
Pfc  T  O  P  'M 
K  0 


